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Abstract 






There are a finite number of identical aircraft each of 
which contain a number of different types of components 
which fail at different rates. In order for an aircraft to 
be operational, all of its components must be operational. 
Each component type has a finite number of spares. Failed 
components are repaired at a single server facility. 
Simulation is used to study the effect of 6 different repair 
scheduling policies. The repair policies are compared on the 
basis of average number of operating aircraft at the end of 
a mission period of one week. It is found that a repair 
policy which first repairs the component of the type with 
the fewest operating components is the best. In particular, 
it is much better than first-in, first-out, and also may 
well improve upon a policy that serves the longest waiting 
line first. A simple spares stockage policy is developed and 
evaluated when the above scheduling policy is in use. 
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I . INTRODUCTION 



This project was defined by and is under the direction 
of Professor Donald P. Gaver of the Department of Operations 
Research, of the Naval Postgraduate School. The objective is 
to study the means by which a maintenance policy might be 
more effective in support of a client. In this study, the 
client is a detachment of U.S. Navy P-3 aircraft. 

This thesis presents a comparison of six different types 
of repair policies which a server may implement on a queuing 
line of different components. These components are 
subsystems of larger, identical systems (eg. radios and 
aircraft). Each of these components are mission-essential to 
the aircraft and fail independently of each other. There are 
nt spares allocated to each component i, i = 1,...,I. The 
repair facility and the spare provisioning constitute the 
maintenance effort of the detachment. The repair facility 
has a single server who must decide on priorities for 
repairing the different types of failed components 
(Figure 1 ) . 

Simulation is used to study six repair policies for the 
server. The first priority policy is First-In, First-Out; 
the second policy serves the most numerous component type in 
queue first; the third serves first the largest product of 
the numbers of component types and individual component 
traffic intensities; the fourth's criteria is to serve the 
least frequent failure component first (smallest failure 
rate); the fifth repairs the most frequent failure component 
first (largest failure rate); and the sixth scans the 
current operating inventory for each component and repairs 
the component with the fewest operating. The six policies 
are compared on the basis of the average number of aircraft 
that have all I mission critical components up after a 
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mission of length T. The results of the simulations are 
presented in tabular form in Appendix B. 
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Figure 1.1. System Diagram with Repair Facility 



Maintenance and Material Management data from the Naval 
Aviation Logistics Data Agency (3-M/NALDA) are used for 
estimating failure rates and repair rates. In addition, 
extensive interviews with personnel at the Aviation 
Intermediate Maintenance Department at Moffett Field, 
California; and Naval Air Systems Command at Washington 
D.C., supplemented the data. The interviews indicated a need 
to look at the small detachment (less than six aircraft) 
problem as well as the short detachment problem (less than 
90 days ) . 

The Aviation Detachment consists of a small number of 
aircraft, each of which contains large numbers of 
components, sub-components , sub-sub-components , and so forth. 
Each of these systems contribute to the availability of a 
given aircraft's readiness for flight. Each aircraft, in 
turn, contributes to the detachment's ability to carry out 
its mission for some duration, which is typically not a long 
period of time. 
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Simplifying assumptions are made in order to reduce the 
size and complexity of the analysis, aid in model 
verification, and provide understandable results. The only 
equipment considered for this study were pieces of 
electronics gear which were designated as mission-critical 
and failure-prone. The length of a mission is measured in 
flight hours and approximates one week's worth of flying. 
Resupply of the detachment occurs at that time; all parts 
are replenished. 

The simulation models were implemented using FORTRAN 77 
programs that are listed in Appendix A. The user inputs 
failure rates of specified parts, total mission time before 
resupply, total components, and total initial spares per 



component . 


The 


output displays 


the average number 


of 


"up 


aircraft" 


at 


the time of 


mission completion, as 


well 


as 


distr ibuti 


onal 


information . 


A detailed discus 


sion 


is 



provided in Chapter III. 
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II. NATURE OF THE PROBLEM 



The detachment of small numbers of P-3 Orion aircraft 
by the U.S. Navy to remote sites around the world is a 
current method of operation by Fleet Commanders. These 
detachments are of relatively short duration, about six to 
eight weeks. They are required to be self-sustained or 
sustainable through air resupply. Small and measured amounts 
of replacement components arrive by P-3 or C-130 aircraft. 
The ability of these detachments to perform their mission is 
heavily dependent on (a) spare parts flown in initially and 
then at specified intervals, and (b) maintenance repair 
policy and capability. 

In the case of the first, (a), considerable effort has 
been expended to address the spare part requirements for P-3 
detachments to remote areas. However, the bulk of these 
consider larger detachments than three aircraft (Ref 1], and 
longer duration than six weeks. The case of provisioning 
smaller, shorter duration groups is most frequently achieved 
by senior maintenance personnel using corporate knowledge 
and experience to derive the requirements. In the case of 
the second (b), little information is available on the study 
of the effects of different repair policies. General 
experience by the author indicates that repair policy is a 
function of the current maintenance administration (policy 
as it relates to which part to fix when). 

The overall goal of the Chief of Naval Operations is to 
achieve at least seventy-two percent of fully mission 
capable aircraft in a squadron. In whole numbers of 
aircraft, this translates to two of three in a small 
detachment. But this goal belies the fact that for a small 
group of planes to accomplish its mission, clearly all the 
aircraft must be available most of the time. 
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The equipment selected for the present analysis was 
avionic gear which had, (a) mean time between failures 
short enough to ensure a reasonable chance of failure during 
a period of detachment, (b) was determined to be mission 
essential by proper authorities, and (c) was repairable or 
replaceable by the detached personnel. Each piece of gear 
was considered to fail independently of the others. The 
overall mission availability of an aircraft is modelled as a 
series system. If all the components are up, then the 
aircraft is up. 



Availability 





Figure 1.2. Minimum Path Representation of N 
Component System 



The failure interarrival times at the queue of all 
components are assumed to be exponentially distributed. The 
service times are also assumed to be exponential. For the 
first series of simulation runs, the failure and service 
rates were the same for all parts. For the latter simulation 
runs, individual failure and repair rates for components are 
used . 
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III. MODEL DESCRIPTIONS 



A. FAILURE MODELS 

There are a variety of approaches to the modelling of 
component supply in a multi-component system with spares 
allocation. Generally speaking, as components fail and 
spares are utilized for replacements, spare part inventories 
can be entirely depleted. Any additional failures which 
occur beyond that inventory can result in systems as a whole 
being unavailable as they await maintenance. The rate at 
which parts arrive at the service center is a constant until 
that point is reached because the number of parts in use at 
any given moment, based on the operating aircraft, is 
constant. When the number of operating systems (aircraft) 
starts to drop, so do the corresponding arrival rates. A 
model incorporating this effect is referred to as the 
decreasing arrival rate model. This is a simple situation to 
simulate, and to obtain results for. However, it is not easy 
to derive analytical results for it. A more analytically 
tractable model is to assume arrivals of a part type to the 
service center form a Poisson process. This will be 
hereafter referred to as the constant arrival rate model. 
This latter approach is rationalized as follows; the 
remaining systems are required to increase their work load 
(flying hours) to compensate for the loss of a system(s), 
ergo more strain on the remaining systems. In addition to 
this, though, a higher (constant) arrival rate would 
represent a conservative approach to calculating spare 
requirements. A higher arrival rate would yield greater 
spares allocation in provisioning (most provisioning models 
use failure rates as the prime method of specifying spare 
requirements ) . 



6 



Initially, two simulation models were written to compare 
the two models for numbers of failures: constant, and 
alternatively decreasing, arrival rates. Ideally, if the 
constant rate simulation results did not depart excessively 
from the actual decreasing rate simulation, then it would be 
desirable to use the constant rate model because it could be 
studied analytically. 

Monte Carlo simulations were written in Fortran 77 and 
utilized a proven random number generator (LLRANDOM II). The 
simulations modeled aircraft components failing and being 
replaced by available spares; failed components were not 
serviced by a repair facility. The simulations used event 
stepping from failure to failure with exponential times 
between failures. To facilitate the analytical analysis, the 
number of aircraft was set at three, the total number of 
types of components per plane at two, and the number of 
spares for each type of component at one. The failure rate 
of component 1 is 0.02 failures per hour; the failure rate 
of component 2 is 0.0143 failures per hour. Each aircraft 
requires both component types to be in operation so that it, 
in turn, may operate. The generated output was the average 
number of up aircraft at the end of a specified time. The 
constant arrival model was verified through calculation 
(Appendix A). This simulation was then modified slightly to 
cause the failure rate to decrease when whole systems 
dropped off line to create the second model. Both 
simulations had five hundred replications. The results are 
listed below in Table III-l. 



TABLE III-l. RESULTS OF MONTE CARLO SIMULATIONS OF 

DECREASING, CONSTANT, AND ANALYTICAL 
FAILURE RATE MODELS 

Expected 

Initial Initial Time in Up 

Model Aircraft Spares Hours Aircraft 



Decreasing 


3 


1 


120 


1.022 


Constant 


3 


1 


120 


0.022 


Analytical 
( Constant ) 


3 


1 


120 


0.0201 
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The conclusion drawn from the simulations is that the 
constant failure model does not describe well the situation 
of the failure rate for the components decreasing as whole 
systems drop off line. However, when whole systems or 
aircraft drop off line for maintenance, the remaining 
aircraft must fly the same work load. This results in 
subjecting the components on the remaining planes to more 
hours in operation. Therefore, the assumption that demand 
rates for spares decrease as whole systems drop off line is 
not entirely correct either. This situation of increased 
hours of operation on remaining aircraft is not easily 
quantifiable and might be consequential. The disparity 
between the models may be less dramatic for this reason. For 
this reason, the constant model is more commonly used in 
application (e.g., in various METRIC models at the RAND 
Corporation) [Ref. 2). Reality is probably a middle ground 
between the two. However, the constant arrival rate model 
will be used in this with the caveat that it produces a 
conservative approach to provisioning. 

B. QUEUING MODELS 

The repair scheduling policy that a maintenance facility 
uses directly affects the available Inventory of spares. How 
dramatically a given policy may influence systems as a whole 
is of concern to any maintenance supervisor. Six repair 
simulation models were written in Fortran 77 to study the 
effects of different repair policies on keeping whole units 
or systems operating. Model algorithms are listed in 
Appendix C. Example Fortran code is listed in Appendix D. 

1. Simulation Background 

The baseline simulation models an aircraft 
detachment located at a remote site. The detachment consists 
of several aircraft, the collection of which are serviced by 
a single maintenance facility. This facility has a single 
server who can work on a single part at a time. The server 
draws his parts from a single queue. As soon as he completes 
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a repair, he immediately commences work on the next part if 
a backlog exists. Each aircraft contains several component 
types, each of which has a spare allocation (although the 
allocation may be zero). Each aircraft requires one of each 
component to operate. The detachment is expected to conduct 
operations for a period of time, and to be self-sufficient. 
A resupply occurs at the end of this period and all spares 
are replenished. The primary measure of effectiveness is 
average number of operating aircraft at the time of 
resupply; variance of the number of operating aircraft is 
also tabulated; the empirical distribution could be 
tabulated if desired. For this study, resupply occurs at the 
end of one week or approximately one hundred and twenty 
flight hours. 

The assumptions which were incorporated in the 
simulation are as follows: the interarrival times of 
component failures are described by an exponential 
distribution; each component has a unique failure rate and 
fails independently of other components; failed components 
are instantaneously removed from the aircraft, replaced with 
available spares, and placed in either the repair queue or 
directly with the repair server; the arrival of failed 
components to the repair facility is described by a Poisson 
process with a rate equal to the sum of all component 
failure rates times the number of detached aircraft; 
cannibalization or the interchange of components to maximize 
working whole aircraft will occur; an aircraft needs all of 
its components to operate and anything less constitutes a 
down aircraft. Consequently, arrivals of failed components 
to the repair facility form a Poisson process with constant 
rate which does not change if the number of whole systems 
change. However, if no aircraft are operating, then the 
arrival rate is 0. 
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is of type i with 



A component which arrives for repair 
probability: 

Ai / (Ai + As + A» + * * * + Ax ) 

Component i's service time is described by an exponential 
distribution with parameter Hi (mu); no balking occurs at 
the queue; the queue has potential length equal to the total 
number of components in the system; and the repair facility 
is capable of repairing all components. 

2 . Simulation Structure 

All simulations evolve from a basic "failure-repair" 
model with a First-In, First-Out service policy. These 
repair models simulate continuous time systems by using 
"event stepping" of a simulation clock time from failure to 
repair, failure to failure, repair to failure, repair to 
repair events. 

All simulations are structured with a main program to 
read input files, generate failures, control simulation 
clock and run parameters, and print statistics; a queue 
subroutine for repairs and policy decision criteria; a 
statistical subroutine to assimilate run data and generate 
output statistics; a random number subroutine to call the 
IMSL library and the LLRANDOM II package to generate arrays 
of pseudorandom uniform variates; and, in all but the 
baseline simulation (First-In, First-Out policy), a priority 
subroutine to assign a component priority based on decision 
cr iter ia . 

Input parameters are total mission duration (in-flight 
hours), total number of components per aircraft (parts), 
total number of spares allocated per component, and total 
number of replications of the simulation. The input file 
contains specific failure and repair rates for each 
component . 

The output from a simulation run contains the arrival 
rate (lambda) at the queue; service rate and traffic 
intensity in those cases where service rates for all parts 
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were set equal; expected numbers of operating aircraft at 
mission termination, variance and standard deviation; a 
tally of the number of Up Aircraft at the end of each run, 
for all runs; average maximum queue size; and average wait, 
in hours, by a part in the queue. 

The repair policies of the various simulations differ at 
the waiting line of the queue. Policies adjust part 
positions in line for service by current repair decision 
cr i ter ia . 

3 . Simulation Techniques 

Pseudorandom number generation is accomplished by 
utilizing a proven pseudorandom number package resident in 
the computer library titled LLRANDOM II; it was developed by 
Dr. P.A.W. Lewis of the Naval Postgraduate School (to check 
for the effect of the starting seed) (Ref. 31. A subroutine 
calls LLRANDOM II three separate times, with three separate 
seeds, to generate three distinct arrays to store the 
uniform random variates. These uniform random variates are 
used to generate failure interarrival times, repair times, 
and determine which part has failed. 

Failure interarrival times and repair times are 
exponential and are calculated using an inverse transform 
method. By incorporating this method into the simulations, 
any of a variety of distributions could be selected for the 
repair times (using the memoryless property of the 
exponential distribution heavily in generating failure 
times). The algorithm states that given a cumulative 
distribution function, F(x) = P(XSx). F ( X ) is uniformly 
distributed over the interval zero to one. By equating a 
uniform random variate to the CDF, F(X), and solving for the 
inverse of F(X), random variates with the distribution of 
F ( X ) can be generated. (Ref. 4] 
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Sample data on the number of up aircraft at mission 
termination is assimilated in the STATS subroutine. Sample 
statistics are then calculated where: 

Xu = Number of up aircraft at mission's end on run j 
n = Total number of runs 

Qi = Maximum queue size during simulation run i 
Wi = Total amount of time (hrs) spent waiting for 
maintenance by the i th arrival at the queue 
m = Total number of parts which waited in queue for all 
runs 

n 

X = Etnumber of up aircraft} = E Xu / n 

j=l 

Var { X } = Variance of X = E (Xt - X) 2 / (n - 1) 

i = 1 

S.D. = Standard Deviation of X = (Var X)* 9 

m 

Wq = Average wait in the queue = E^Wi / m 



n 

Q = Average maximum queue size = E^Qi / n 

4 . Simulation Walkthrough 

A simulation run begins with the initialization of 
mission parameters such as stock levels of components 
(components in use plus spares), mission clock, subroutine 
parameters for the queue, various flags and counters. The 
initial failure interarrival time is calculated and compared 
to an initial repair time (initially set high). The lesser 
of the two times determines the next event and the mission 
clock is advanced to that exact moment. Stock levels are 
adjusted when the events occur; repairs increment, failures 
decrement stocks. In the case of failures, after the 
appropriate stock is reduced, a component is instantly 
placed in the repair facility where it may either be served 
immediately or join a queue (repair policy will determine 
what position in line the part assumes). After a failed part 
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has arrived at the repair facility, all component stock 
levels are scanned to ensure at least one aircraft remains 
in operation. If all aircraft are down, failures cease, the 
mission clock is advanced to the next repair time, and the 
repair event occurs. If at least one aircraft is up, failure 
and repair events continue as before. In the case of a 
repair event, after the appropriate stock is incremented, 
the next part in line is immediately brought into service 
and a new repair time is calculated. The component stocks 
are scanned again to ensure at least one aircraft is up. If 
the event where previously no aircraft were up but because 
the correct type of part was fixed one now exists, a failure 
interarrival time is then calculated and event comparison 
transpires as before. If the all aircraft are still down 
after a repair, the clock again advances to the next repair 
time and that event occurs. 

The mission clock is allowed to be advanced from event 
to event until the scheduled time exceeds the predetermined 
mission duration. When this value is met or exceeded, the 
simulation run stops, data is gathered on those left waiting 
for service in the queue (waiting time) and the number of 
operating systems is recorded. The run number is compared to 
the total number of replications value input at the start. 
If more runs are required, a mission profile is 
reinitialized and the process repeats itself. Otherwise, 
statistics described in previous sections are computed and 
then displayed (at the terminal). 

5 . Flrat-In, Firat-Out Model 

This model is the basic simulation from which all 
others were extended. The First-In, First-Out or FIFO is a 
common form of maintenance scheduling policy. When a part 
arrives in the queue it is placed at the end of the line. 
When the server becomes available to work, the part at the 
head of the line or first to arrive, is served next. There 
are never any interruptions once service starts. 
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6 . Dynamic 1 Model 

This policy adopts the following criteria: use a 
First-In-First-Out routine until such time as there are two 
or more of a given type of part in line, at which time the 
more numerous part type receives first priority for repair. 
For example, there are three types of parts in the queue 
awaiting maintenance, parts A, B, and C. There are two part 
A's, three part B's, and one part C. Parts B would move to 
the front of the line and be serviced first as long as they 
remain the most numerous. 

7 . Dynamic 2 Model 

This service rule tells the repairman to count the 
number of parts of each type of component and multiply this 
by the individual traffic intensity of the component, then 
pick the largest value. The rule in the form of equations is 
as follows: 

p t = Xt / Mi = Traffic intensity of component i 

Priority = Maxl^i x Ni(t)> 

Where : 

Ni(t) = the number of parts of type i in the service 
line at time t 

Hi = the service rate of component i 
/i = the failure rate of component i 

8. Dynamic 3 Model 

This scheduling policy determines service priority 
based on the current available stock of operating parts for 
each type component. After the repairman scans the stock 
levels, he reorders the parts in the service line to favor 
the one with the lowest operational inventory. 

Failure Rate Priority-Low 

This model assigns priority of service based on the 
failure rate of the type part. Those with the smallest 
failure rate, or least frequently failing, are given first 
service regardless of their arrival time. When the server 
scans the queue for the next part to fix, he compares the 
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individual part lambdas to select the next repair. If 
multiple parts are of the same type, a FIFO policy is used. 
The abbreviation FRP-Low is used throughout the text. 

10 . Failure Rate Pr lor i tv-Hiqh 

This model is similar to FRP-Low with the following 
modification; priority is assigned to the most frequently 
failing or highest failure rate part. If multiple numbers of 
the same part are being considered, a FIFO policy is used. 
The abbreviation FRP-High is used throughout the text. 

C . DATA 

The primary source for failure and repair rate data in 
this study is the NALDA (Naval Aviation Logistics Data 
Agency) data base. The reason this source is used is because 
it employs extensive error checking algorithms in its 
database. It is considered the best source of clean data 
(error reduced) for P-3's by the experts at Naval Air 
Systems Command and in the fleet. The information it 
supplied was supplemented by personal interviews by the 
author with maintenance supervisors at the AIMD, N.A.S. 
Moffett Field, California. 

The interviews cautioned that even though extensive 
error checking algorithms are employed by the data base, 
they were not error free. Precise values for the Mean Time 
Between Failures (MTBF) and the Mean Time To Repair (MTTR) 
would require extensive data analysis. Reasonable 
approximations, however, could be easily extracted. 
Interviews provided an understanding of the sources of the 
data entries and the means to extract reasonable approximate 
values . 

The focus of this thesis is repair policy evaluation. 
The author sought out values for failure and repair rates 
which would test repair decision criteria with a range of 
values and be representative of mission essential avionics. 
The final values are listed in Table III-2. The complexity 
of the P-3 avionics suite makes it necessary to approximate 
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failure and repair rates for components. This is done 
because a component is composed of subcomponents, which are 
in turn composed of sub-subcomponents. Data information is 
gathered on the subcomponent level but requirements are 
derived on the whole component level. A separate study can 
be done to determine the component failure and repair rates 



as a whole. This thesis 


required 


only approximate 


values . 


Failure rates 


are derived 


from the 


mean 


number of flight 


hours between 


failures (MFHBF). The 


repair 


rates are 


derived 


from dedicated 


maintenance man hours (DMMH) 


at the 


AIMD, or 


intermediate maintenance 


level. It 


is important to note that 


the precise numerical values are not important as 


long as 


the relative 


magnitudes 


are correct and are reasonably 


representative 


of actual 


values. 








TABLE III- 


2. FAILURE 


AND REPAIR 


RATES 






CortiDonenti 


Lambda i 


MFHBF 


Mui 


DMMH 




1 


0.0200 


50.0 


1.000 


1.0 




2 


0.0143 


70.0 


0.417 


2.4 




3 


0.0180 


55.5 


1.000 


1.0 




4 


0.0026 


389.0 


0.333 


3.0 




5 


0.0193 


51.8 


0.666 


1.5 




6 


0 . 0219 


145.8 


1.333 


0.8 




7 


0.0010 


1000.0 


0.200 


5.0 




8 


0.0110 


90.9 


0.400 


2.5 




9 


0.0300 


33 . 3 


2.000 


0.5 




10 


0.0060 


166 . 7 


0.444 


2.3 
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IV. SIMULATION RVM£ 



A. BACKGROUND 

The six repair policies were exercised through a series 
of different case studies in which input parameters were 
varied identically for all six policies. The Measure of 
Effectiveness (MOE), average number of aircraft in operation 
at the termination of the mission for each repair policy, 
was then compared. Due to the experimental nature of these 
repair policies, a control environment was needed to provide 
a benchmark of performance. Case studies 4 and 5 were an 
attempt to achieve this effect. Their model inputs were set 
equal across all component types (a more detailed discussion 
follows). Embellishments of the system operating environment 
were conducted in Case 6 and evolved further in Cases 1 
through 3. These refinements subjected the models to ever 
more variable and complex system parameters. 

The input parameters which were altered were: numbers of 
aircraft (planes); spares allocation for each component 
type; and individual component failure and repair rates. The 
total number of component types contained on each aircraft 
was held constant at 9 for all repair policies and cases. 
All simulations had five hundred replications. 

1 . Spares Allocation 

Two different schemes for changing spare part 
stockage levels were devised. The first scheme is referred 
to as the "Spare to Plane Ratio" method. Initial spare 
provisions are set equal for all component types (i.e. 
component A’s spares are initialized at 12 parts, as are 
component B’s, as are C's, etc). As the name of the method 
indicates, the amount which the spare levels are varied is 
in proportion to the number of aircraft in the detachment. 
This ratio changes from 2, to 1.5, to 1, to 0.5, to 0.1666 
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if there are six detachment aircraft 



(1/6). For example, 

(Planes), all components spares are varied from 12 (spare to 
plane ratio of 2), to 9 (1.5), to 6 (1), to 3 (0.5), to 1 

(0.166). These provisioning changes constitute the 
differences between Cases 4a-4e, 5a-5e, and 6a-6e. 

The second scheme is referred to as the "K-Standard 
Deviation" method. Predictions of the average number of 
failures a component will experience during a mission period 
is a measure of the demand for spares. The Poisson process 
is employed to estimate the mean number of spares each 
component would require during a mission period. The 
sensitivity of this estimate of the demand for spares is 
tested by adding to or subtracting from the mean K standard 
deviations of spares (from the appropriate Poisson 
distribution). For example, if the mean number of spares 
plus two standard deviations were stocked, the demand for 
parts could be expected to be met 95% of the time (without 
repair). The provisioning calculations proceed as follows: 
Expected Demand for Part i = (a x t x At ) = ElXil 
Standard Deviation of Demand = ( ElXil )•* 

Spare Policy Component i = ElXil + K x ( ElXil )•* 

Where : 

a = Number of aircraft 
t = Total mission length 
X i = Failure rate of Component i 
K = Standard Deviation factor 

The result of this policy for components having different 
failure rates is stock levels which are different for each 
component type. The most frequently failing component 
receives the most spares. Variations of the value of K 
determine the difference in Cases la-lh, 2a-2h, and 3a-3h. 

2 . Failure and Repair Rates 

Two separate parameter schedules are used in the 
experimentation process. In the first, failure and repair 
rates are set equal for all component types (i.e.. Component 



18 



A's failure rate is 0.02 failures per hour and its' repair 
rate is 0.84 units per hour, as are Component B's, etc). 
This input schedule is used for Cases 4a-4e, and 5a-5e. The 
second schedule uses individual failure and repair rates 
listed in Table III-2 (Failure rates = Lambda!, Repair rates 
= Mui). This parameter schedule is used for Cases la-lh, 2a- 
2h, 3a-3h, and 6a-6h. 

3 , Output 

The output of each case is given in Appendix B. The 
listed results include the Mean number of operating aircraft 
at mission termination (Average Up Planes), the Variance of 
the number of up aircraft, the Standard Deviation of the 
number of up aircraft, the Standard Error of the Mean, and 
the fraction of Up Planes operating at mission termination. 

The number of operating aircraft at the end of a given 
run is determined by analyzing the number of operating parts 
of component type i, i = 1,...,I. The component with the 

minimum number of operating parts determines the maximum 
number of operating planes. The Mean number of operating 
aircraft is calculated from the numerical average of Up 
Aircraft at the end of all replications (runs). This is the 
common sample mean: 

500 

Mean Up Aircraft = E Xi / 500 = X" 

i = 1 

Where : 

Xi = Number of Up Aircraft on run i 
The Variance of the run sample and the Standard Deviation 
are : 

500 

Variance = E ( Xi - X) 2 / 499 = VartX] 

i = l 

Standard Deviation = ( Var[X] )•* 

The Standard Error of the Mean is given by the relation: 

Standard Error = S_ = ( Var(X) / 500 ) • 9 

The Fraction of Up Planes operating at mission termination 
is used as a measure of effectiveness for the policies. It 
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is calculated after each subcase is completed and is derived 
by dividing the Mean number of operating Aircraft ( E [ X ] ) by 
the total number of detachment aircraft. (Ref. 5] 

4 . Random Seeds and Stability 

In order to establish that the simulations had 
reached stability in sampling variance, five of the six 
models were simulated by varying only the random number 
seeds used in generating failures and repairs 1 . Four sets of 
seeds were utilized. Two-way Confidence Intervals for the 
Means (Expected Number of Up Aircraft) were then computed at 
a level of 95%. The distribution of the Sample Means was 
assumed to be Normal in accordance with the Central Limit 
Theorem. The results are listed in Appendix B. and 
demonstrate that the sampling variance had stabilized. For 
any given policy, the Means of any combination of three 
random number sets fell within the confidence interval of 
the remaining fourth set. For any given set of random number 
seeds, the six repair policy Means established a rank order 
by magnitude. This relative ranking between policies did not 
change for all four seed sets. As a result, it was 
determined that reliable output from the simulations could 
be anticipated with five hundred replications. 

B. CASE DESCRIPTIONS 

1 . Case 4-5 

These cases were designed to give a baseline measure of 
how the models would compare in a environment where 
interrelat ional effects between parameters such as failure 
and service rates, and stock levels, could be minimized. 
Case 5 differs from Case 4 in the number of total aircraft. 
Subcases of Case 4 and 5 delineate variations in initial 
spare provisioning as described by the "Spare to Plane 
Ratio" method. Inherent to the background of both cases are 



^-The sixth model. Dynamic 3 was added after this 
exercise. Stability of variance for it was assumed based of 
the other models' results. 
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failure and repair rates, and total numbers of components 
per plane (9). Failure and repair rates are fixed at equal 
values (0.02 and 0.84 respectively) for all types of 
components. These rates represent average values for the 
listing in Table 1 1 1 — 2 . 

2 . Cases 1-3 

In these cases, the failure and service rate for the 
components differ. These cases were designed to analyze how 
sensitive the various models were to a dynamic environment 
where several of the parameters are varied (and may have 
potential interrelated effects). For example, high failure 
rates and high service rates in conjunction with large 
quantities of spares may result in the server dedicating a 
significant portion of his effort towards a single type 
part. As before. Cases 1, 2, and 3 differ in numbers of 

detached aircraft. Case 1 contains three; Case 2 contains 
six; Case 3 contains twelve. Subcases delineated spare 
provisioning as described by the "K-Standard Deviation" 
method. Cases 1, 2, and 3 differ in their selection of 

values of K. The range of K is +1 to -4. Negative values 

of K are selected to explore the effects of decreasing spare 
levels. The actual spare allocations as a function of case 
number are listed in Table B-3 of Appendix B. An excerpt 
from that Table is given below as an example. 

TABLE IV-1 MODEL PARAMETERS: SPARE ALLOCATION OF 

COMPONENTS FOR SIMULATION CASES 

Casa la Case lb Case lc Case Id 



K-Factor 


+ 1 


0.0 


-0.5 


-1 


PART 1 


10 


7 


6 


5 


PART 2 


7 


5 


4 


3 


PART 3 


9 


6 


5 


4 


PART 4 


2 


1 


0 


0 


PART 5 


10 


7 


6 


4 


PART 6 


11 


8 


6 


4 


PART 7 


1 


0 


0 


0 


PART 8 


6 


4 


3 


2 


PART 9 


14 


11 


9 


8 
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Unique component failure and service rate are input from 
Table 1 1 1 — 2 . Total components per aircraft are nine. 

3 . Case 6 

This was an intermediate step from Cases 4-5 and Cases 
1-3 for parameter selection. Spare stockage levels employed 
the "Spare to Plane Ratio" method as per Cases 4 and 5. This 
determined the variation in the subcases (6a-6e). Unique 
failure and repair rates are initialized at the values 
listed in Table III-2, and as per Cases 1-3. The number of 
component types per aircraft is unchanged at nine. 
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V. RESULTS OF CASE STUDIES 



The purpose of the case studies is to compare different 
repair policies using as the primary MOE average number of 
aircraft up at the end of a mission of finite length. A 
repair policy substantially more responsive to shortages 
then a simple First-In, First-Out routine should stand out. 
The Dynamic 3 policy consistently achieved this result. 
Conversely, the Failure rate priorities are consistently 
poor performers (relative to the FIFO benchmark). The other 
two Dynamic policies are sensitive to the spare provisioning 
policy. Recall the fact that both Dynamic 1 and Dynamic 2 
set priority for service based on total numbers of component 
types awaiting maintenance. High-failure rate parts, which 
are also provisioned more heavily, will tend to consume the 
available server's capacity. Those parts which are stocked 
at lower levels (based on lower failure rates) are relegated 
to the "back of the line". It is the effect of all the lower 
rate failure items being placed at the end of the line which 
causes the Dynamic 1 and 2 to drop in performance (when 
there are unequal numbers of spares for each component). 
There is a cross-over point for these two models relative to 
the FIFO policy, that point being when the spare provisions 
are depleted to the point where stock levels begin to 
equalize (K factor i -2.25). More detailed case by case 
discussions are provided in this chapter. 

Case study results are listed in Table B-l, Appendix B. 
The graphical summaries contained in this chapter provide a 
visual means to summarize the effectiveness of each of the 
repair policies relative to each other. Since Cases 4-5 
provide a baseline measure for the models, they will be 
discussed first, followed by Case 6, and then Cases 1-3. 
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A. CASES 4-5: EQUAL FAILURE AND REPAIR RATES 

The equality of failure rates for each component and 
repair rates for each component should demonstrate the 
relative capability of each repair policy to replenish spare 
stock levels and keep airplanes operating. A significant 
difference between policies can be determined by comparing 
the Mean and the Standard Error for each data point from 
Appendix B and calculating Normal confidence intervals for 
the mean (also by recalling that for Normal distributions, 
95% of a population distribution should lie roughly within 2 
standard deviations of the mean). By analyzing the data and 
graphs (Figures 5-1, 5-2), for Cases 4-5, the Dynamic 
policies (1, 2, 3) are consistently more able to keep 
aircraft operating (stocks replenished). 
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Figure 5.1. Fraction of Operating Planes as a Function 
of Spares to Planes Ratio (Case 4). 

The difference between them and the FIFO routine is 

significant at the point where the spare to plane ratio is 

equal to 1.5. Since failure and service rates are equal for 



all components. Dynamic and Dynamic 2 are in essence 
following the same policy as Dynamic 3; that being, repair 
the lowest stock level first (or conversely the most 
numerous in the queue). 




Figure 5.2. Fraction of Operating Planes as a Function 
of Spares to Planes Ratio (Case 5). 

The actual differences in the policies are attributed to 

variation within the simulation. The means for the Dynamic 

policies were within 1.96 standard deviations of each other 

(95 % Confidence Interval). The Failure rate priority 

policies were not considered in Case 5 due to their poor 

performance in Case 4. In Case 4, they demonstrated 

performance which was as much as twenty five percent less 

effective than FIFO. The difference between the Failure rate 

priority routine and FIFO is because FRP-L and FRP-H assign 

priorities one through nine to each of the parts. Since all 

the rates were identical, the priority assigned was one for 

component one, two for component two, etc. Also, since the 
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rates were identical, there is no difference in component 
priorities in FRP-L and FRP-H. 

B. CASE 6: DIFFERENT FAILURE AND REPAIR RATES; IDENTICAL 

SPARES STOCKING 

These simulations address the issue of how the policies 
respond to variable failure and service rates with a 
constant part stockage. Individual component failure and 
repair rates are given by Table 1 1 1 — 2 . Again, since the 
numbers of spares is the same for all parts. Dynamic and 
Dynamic 2 still reflect the Dynamic 3 policy (numbers in the 
queue relate to the stock levels which relates to Up 
Planes). The Dynamic policies are still consistently better 
than either FIFO or the Failure Rate policies (Figure 5-3). 
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Figure 5.3. Fraction of Operating Planes as a Function 
of Spares to Planes Ratio (Case 6). 

Again, the point of significant departure for the Dynamic 

policies and FIFO occurred at a plane to spare ratio of 1.5. 

The Dynamic policies' means fell outside of a 95% confidence 

interval for the FIFO policy. The differences between the 
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Dynamic policies is attributed to variation in the Mean 
(within a 95% confidence interval). The Failure rate 
policies showed improvement. These two policies are now 
dissimilar because failure and repair rates are different. 
Between the two, FRP-H performs significantly better than 
FRP-L (based on a 95% confidence interval about the means). 

C. CASE 1-3: DIFFERENT RATES, DIFFERENT STOCKS 

In these cases, we see a change occur in the abilities 
of Dynamic 1 and Dynamic 2 to replenish stock levels. 
Initial stock levels are no longer equal for all parts. They 
vary for each component based on failure rates. High failure 
rate components are supplied the most spares; low failure 
rate components the least. Since Dynamic 1 and Dynamic 2 
repair the component having the most numerous failed parts 
first, the issue of replacing the most needed operating 
stock is skirted. High-failure items can fill the queue and 
still maintain sufficient operating units on all aircraft. 
Other units may fail; however, if they fail in insufficient 
numbers, they may never be served. For example, in a case in 
which the expected demand (failures) is stocked for 
Component 1 and 8 (K = 0 for 6 planes), their respective 
spares are 14 and 8. Component l's failure rate is 
approximately twice eight's (0.02 vice 0.011). Therefore, 
Component l's arrive at the repair facility twice as often. 
If there are 14 Component Is waiting for repair and 13 
Component 8s waiting for repair, then the Dynamic 1 and 2 
policies would choose Component 1 to repair next even though 
there are still 6 Component Is operating and only 1 
Component 8 operating. This policy would result in only 1 
plane operating. A comparison of Case 6 with Case 2 will 
demonstrate what happens when the spare provisioning policy 
is changed from equal spares for all components (Spare to 
Plane Ratio method) to a failure-based system (K-Standard 
Deviation method). The Dynamic 3 policy is concerned only 
with those components which are in most demand (lowest 
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operating stock level), regardless of the number awaiting 
service. It doesn't waste time on repairing those failed 
parts that have sufficient spare stock levels at a given 
moment. Dynamic 3 was consistently the best of all policies 
(Figure 5-4 , 5-5, 5-6 ) . 




Figure 5.4. Fraction of Operating Planes as a Function 
of K Standard Deviations of Parts (Case 1). 

Once again, the Failure policies do not to perform well. In 

Case 1 with 3 planes, there is a tight bunching of all the 

policies. This is due to the fact that the server is not 

very busy with 3 planes and as a result, the form of the 

repair policy is not critical. As the number of planes 

increase, the policies diverge; the server becomes less able 

to rejuvenate stocks in the face of higher demands. 
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Fraction of Operating Planes as a Function 
of K Standard Deviations of Parts (Case 2). 
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of K Standard Deviations of Parts (Case 3). 
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VI, RESULTS AND CONCLUSIONS 



This thesis has presented a comparison of six repair 
scheduling policies. The Dynamic 3 repair model, which 
schedules service in favor of the component with the lowest 
operating stock, yields a significant improvement of 
available aircraft at the end of a resupply period when 
compared to a simple First-In, First-Out policy. A key 
factor in repair policy scheduling is the current inventory 
of available operating components. The other models which 
were studied in this thesis failed to address this issue 
directly. As a result, when spares were stocked in unequal 
numbers (the most likely real world situation), the other 
policies (Dynamic, Dynamic 2, FIFO, FRP-L, FRP-H) did not 
perform as well as Dynamic 3. 

The overall goal of the Chief of Naval Operations is to 
achieve at least seventy-two percent of fully mission 
capable aircraft. The maintenance support of a squadron and 
the supporting airwing are responsible for accomplishing 
that goal. It is clear from the studies conducted in this 
thesis that significant differences in aircraft availability 
can result simply by the manner in which a maintenance 
effort schedules its repairs. A policy such as Dynamic 3 can 
assist in achieving the CNO goal better than FIFO. Spare 
stock levels are maintained more effectively regardless of 
the provisioning policy or the workload in the repair 
facility. However, a simple and sensible stocking policy 
increases effectiveness even more; the K-policy illustrates 
this fact. Future study into other repair and stockage 
policies is recommended to study their effects on system 
availability. Other repair policies may look into the 
effects of different distributional assumptions for failure 
and repair times, as well as repair scheduling. It is also 
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necessary to consider the effect of breakdowns of the repair 
facility itself: How can these be best accommodated? 
Finally, it is desirable to reduce the computational effort 
needed to evaluate repair and stockage policies either by 
improving simulation efficiency by Monte Carlo "swindle", or 
by their replacement by analytic approximations. These steps 
are under current examination. 
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APPENDIX A 



A. THE PROBABILITY CALCULATIONS OF UP AIRCRAFT 

The following calculations are used to verify the 
constant arrival rate simulation model with no repair 
f ac i 1 i ty . 

Parameters : 

Ai(t) = The number of type i components available 
and up during the interval (0,tj 
Fi(t) = The number of type i components which fail 
during the interval (0,t] 

I = The total number of types of components 

a = The total number of aircraft 

n± = The total number of spares of type i component 
Ai = The failure rate for component i 
Each aircraft requires one of each type of component to 
operate. The maximum number of up aircraft at any given 
moment t is equal to the lowest level of up components for 
any given type of part. Therefore, the number of up aircraft 
at any given time t is: 

A(t) = Min{Ai(t), A 2 (t), . . . , A s (t),...,a} (A.l) 

If we use the constant failure rate assumption, 
component i will fail at constant rate equal to a i, even 

if the number of up aircraft is less than a. We assume an 

exponential failure distribution and cannibalization of 
parts. This leads to approximating the failure time of part 
i by a Poisson arrival process. We are interested in the 
probability of a minimum number of aircraft being up at 
time t : 

P {Min # A/C up > k} 



P { A( t ) > k} = it P{ Ai(t)> k} 
i = l 

= P{ # A/C avail. > k} 



(A. 2) 
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since different components are assumed to fail independently 
and based on the number of aircraft hours flown ( 0 , t ] 

The probability that the number of components that fail 
of type i in the interval (0,t] is less than or equal to 
some fixed value j is: 

P{# parts that fail type i in (0,t] < j} = 

j e _a.A ^ (a X ±t ) h 

P{Fi(t)<j} = E — , if j < a+n i (A. 3) 

k = 0 k ! 

P{Fi(t)<j} =1 , if j > a+n i 

Now, recall that the number of failed components plus 
the number of available components equals the total 
components . 

Fi(t) + Ai(t) = a + ni -» A± ( t ) = a + n t - Fi(t) (A. 4) 

We would like to solve for the probability that the number 
of available components of type i at time t is greater than 
or equal to some number k or P{Ai(t) > k}. So substituting 
A. 4. for Ai(t), we solve: 

P (a + ni - Fi(t) > k} -♦ P { Fi ( t ) < a + ni - k} = 

P { Ai ( t ) > k) . 



If we let the number j = a 
previous solution (A. 3) for the 
fail of type i in the interval (0, 

a+ni-k 

P { F i ( t ) < a + n i - k } = E 

L = 0 



+ ni - k and plug into our 
number of components that 
t] : P { Fi ( t ) < j) 

e-^Xi* (aXit) L 



= P { Ai ( t ) > k} 

We now solve for the boundary conditions: 

lower: j < a + n t , j=a + ni-k->k>l 
upper: k < a , since you can not exceed the total number 
of aircraft in the system as an upper bound. 
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In summary, the probability that the number of up 
aircraft is greater than some fixed number k is given by... 



I a + ni-k e-dlAit: (a7ut) L 

P{A(t) > k} = rc E if l<k<a 

i=l L=0 L! 



P{A(t) > 0} = 1 ... the entire probability space 

P{A(t) > a+1} =0 ... you cannot exceed the total number of 

aircraft . 

The expected number of up and operating systems at the 
end of some period t is given by the product of the 
probability of a given value and that value: 

E { A( t ) } = OxP { A( t ) =0 } + lxP{ A( t ) =1} + 2xP { A( t ) =2 } + ... 

The variance of the expected number of up and operating 
systems at the end of period t is given by: 

E{ A( t ) 2 } = OxP { A( t ) =0 } + l 2 xP{ A( t ) =1} + 2 2 xP{ A( t ) } + ... 

Var { A( t ) } = E { A( t ) 2 } - E{ A( t ) } 2 

where: P{A(t)=0} = P{A(t)>0} - P{A(t)>l) 

The constant arrival simulation without repair is 
verified with the parameters: 

I = 2, a = 3, ni = 1, t = 120, x = 1/50, 2 = 1/70 

The probabilities then follow: 

P { A( t ) =1 } = P { Ax ( t ) =1 } x P { As ( t ) > 1 } 

+ P { Ax ( t ) £l } x P { A 2 ( t ) =1 } = 0.01479 
P { A( t ) =2 } = P { Ax ( t ) =2 } x P { A 2 ( t ) > 2 } 

+ P { Ax ( t ) > 2 } x P { A 2 ( t ) =2 } = 0.00273 
P{ A( t ) =3 } = P { Ax ( t ) =3 } x P{ A 2 ( t ) >3} 

+ P { Ax ( t ) > 3 } x P { A 2 ( t ) =3 } = 0.00001 
E { A( t ) } = 1x0.01479 + 2x0.00273 +3x0.00001 
E{ A( t ) } = 0.0201 , E{ A( t ) } 2 = 0.0004 
E{ A( t ) 2 } = 0.0258 
Var { A( t ) } = 0.0254 
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APPENDIX B 



TABLE B 


-1 MODEL 


STATISTICS : 


EXPECTED NUMBER 


OF UP 




AIRCRAFT VS SIMULATION 


CASE NUMBER 






Case la 


Case lb 




Case lc 


Case Id 


Planes 


3 


3 




3 


3 


K-Factor 2 


1 


0 




-0.5 


-1 


Parts/Plane 


9 


9 




9 


9 


Lambda/Mu 


Individual Failure 


, Service rates used 


Dynamic 1 
Mean 


3 . 00 


2.99 




2.93 


2.92 


Var iance 


0.00 


0.01 




0.07 


0.07 


S . Dev. 


0.00 


0.10 




0 . 27 


0.27 


s _ 3 

x 


0 . 00 


0.00 




0.01 


0.01 


% Up A/C 


1 . 00 


0.99 




0.98 


0.97 


Dynamic 2 
Mean 


3.00 


2.98 




2.89 


2.89 


Variance 


0.00 


0.02 




0.09 


0.09 


S . Dev. 


0.00 


0 .14 




0 . 32 


0.32 


S_ 

x 


0 . 00 


0.01 




0.01 


0.01 


% Up A/C 


1.00 


0 .99 




0.96 


0.96 


Dynamic 3 
Mean 


3.00 


2.97 




2.93 


2.93 


Var iance 


0.00 


0.03 




0.06 


0.07 


S . Dev . 


0.00 


0 . 17 




0.25 


0.26 


S 

X 


0.00 


0.01 




0.01 


0.01 


% Up A/C 


1.00 


0.99 




0.98 


0.98 


FIFO 

Mean 


3.00 


2 . 97 




2.93 


2.92 


Var iance 


0.00 


0.03 




0.07 


0.07 


S . Dev. 


0.00 


0.17 




0.27 


0.27 


S_ 

X 


0.00 


0.01 




0.01 


0.01 


% Up A/C 


1.00 


0.99 




0.98 


0.97 


FRP-Low 

Mean 


3.00 


2.96 




2.91 


2.91 


Var iance 


0.00 


0.04 




0.08 


0 . 08 


S . Dev . 


0.00 


0 . 19 




0.29 


0.29 


S_ 

X 


0.00 


0.01 




0.01 


0.01 


% Up A/C 


1.00 


0 . 99 




0.97 


0.97 


FRP_-_Hi_qh 

Mean 


2.99 


2.96 




2.92 


2.92 


Variance 


0.01 


0 .04 




0.09 


0.08 


S . Dev. 


0.10 


0 . 20 




0.29 


0.29 


S_ 

X 


0.00 


0.01 




0.01 


0.01 


% Up A/C 


0.99 


0 . 99 




0.97 


0.97 


2 K-Standard Deviation method 


used. 


see Chapter 


IV 


3 S_ = 

X 


{ ( Var iance/500 ) } • 3 = Standard 


Error 
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TABLE B 



P lanes 

K-Factor 

Parts/Plane 

Lambda 

Mu 

Dynamic 1 
Mean 

Variance 
S . Dev . 



% Up A/C 

Dynamic 2 
Mean 

Variance 
S . Dev . 

S_ 

X 

% Up A/C 

Dynamic 3 
Mean 

Variance 
S . Dev . 

S_ 

X 

% Up A/C 

FIFO 

Mean 

Var iance 
S . Dev . 

S_ 

x 

% Up A/C 

FRP-Low 

Mean 

Var iance 
S . Dev . 

S_ 

X 

% Up A/C 

FRP-Hiah 

Mean 

Variance 
S . Dev . 

S_ 

x 

% Up A/C 



"•Note : 
pract ica 1 ly 



1 MODEL STATISTICS: EXPECTED NUMBER OF UP 

AIRCRAFT VS SIMULATION CASE NUMBER 4 



Case le Case if 



Case lq Case lh 



3 

-1.5 

9 



3 3 3 

-2.0 -2.5 -3.0 

9 9 9 

-Individual Failure rates used 

Individual Service rates used 



2.89 


2.65 


2.03 


1.93 


0.12 


0.32 


0.38 


0 . 39 


0 . 34 


0.57 


0.61 


0.63 


0.02 


0.03 


0.03 


0.03 


0.96 


0 . 88 


0.68 


0.64 


2.87 


2.63 


2.09 


1.89 


0.12 


0 . 33 


0 . 39 


0.31 


0.34 


0 . 56 


0.62 


0 .63 


0.03 


0 . 03 


0.03 


0.02 


0.96 


0 . 88 


0.69 


0.63 


2.91 


2.68 


2.04 


1.91 


0 . 09 


0 . 29 


0 .42 


0.36 


0.31 


0.55 


0.64 


0.59 


0.01 


0.02 


0.03 


0.02 


0.97 


0.89 


0.68 


0.63 


2.89 


2.64 


1.96 


1.92 


0.11 


0 .34 


0 . 38 


0 . 36 


0.33 


0.59 


0.62 


0.60 


0.02 


0 .03 


0.03 


0.03 


0.96 


0.88 


0.65 


0.63 


2.86 


2 .65 


2.02 


1.93 


0.17 


0.37 


0.39 


0.33 


0 . 41 


0.61 


0.63 


0.57 


0.02 


0.03 


0.03 


0.03 


0.95 


0.88 


0.67 


0.64 


2.86 


2.65 


2.05 


1.93 


0.16 


0.35 


0.38 


0.34 


0.39 


0 . 59 


0.62 


0.59 


0.02 


0.03 


0.03 


0.03 


0.95 


0.88 


0.68 


0.64 



None of the policies differ much, either 
or statistically. 
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TABLE B-l MODEL STATISTICS: EXPECTED NUMBER OF UP 

AIRCRAFT VS SIMULATION CASE NUMBER 





Case 2a 


Case 2b 


Case 2c 


Case 2d 


Planes 


6 


6 


6 


6 


K-Factor 


1 


0 


-1 


-1.5 


Parts/Plane 


9 


9 


9 


9 












LaUBDUa 




■ ' inuiviuUai failure 


IT cl l e S US 0u 




Mu 




-Individual Service 


rates used 




Dynamic 1 










Mean 


5.99 


5.93 


5.27 


5.23 


Var iance 


0.01 


0.09 


0.71 


0 . 72 


S . Dev . 


0 . 10 


0 .32 


0 .84 


0.85 


S_ 

9C 


0.00 


0.01 


0.04 


0.04 


% Up A/C 


0.99 


0.99 


0 . 88 


0.87 


Dynamic 2 










Mean 


5.96 


5.83 


5.03 


5.03 


Var iance 


0.01 


0.25 


1.09 


1.09 


S . Dev . 


0.11 


0.49 


1.05 


1.05 


S_ 

9C 


0.00 


0.02 


0.05 


0.05 


% Up A/C 


0.99 


0.97 


0.84 


0.84 


Dynamic 3 s 










Mean 


6.00 


6 . 00 


5.88 


5.83 


Var iance 


0.00 


0.00 


0.11 


0.19 


S . Dev. 


0.00 


0.00 


0.34 


0.43 


S_ 

3€ 


0.00 


0.00 


0.01 


0.02 


% Up A/C 


1.00 


1.00 


0.98 


0.97 


FIFO 










Mean 


5.99 


5.99 


5.61 


5.53 


Variance 


0.01 


0.02 


0 . 42 


0.65 


S . Dev. 


0.10 


0.15 


0.65 


0 . 80 


S_ 


0.00 


0.01 


0.03 


0.04 


% Up A/C 


0.99 


0.99 


0.93 


0.92 


FRP-Low 










Mean 


5.94 


5.69 


4.97 


4.58 


Variance 


0.22 


1.11 


2.95 


4.09 


S . Dev . 


0.46 


1.05 


1.72 


2.02 


S_ 

9C 


0.02 


0.05 


0.08 


0.09 


% Up A/C 


0.99 


0.95 


0.83 


0.76 


FRP-Hiah 










Mean 


5.96 


5.71 


4.78 


4.69 


Variance 


0.07 


0.66 


2.03 


2.54 


S . Dev . 


0.26 


0.81 


1.42 


1.59 


S_ 

X 


0.01 


0.04 


0.06 


0.07 


% Up A/C 


0.99 


0.95 


0.80 


0.78 



“Note: Dynamic 3 looks very good in Cases 2c and 2d. 
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TABLE B-l 



MODEL STATISTICS: EXPECTED NUMBER OF UP 

AIRCRAFT VS SIMULATION CASE NUMBER 





Case 2e 


Case 2f 


Case 2a 


Case 2h 


Planes 


6 


6 


6 


6 


K-Factor 


-2 


-2.5 


-3.0 


-3.5 


Parts/Plane 


9 


9 


9 


9 


Lambda 












I na lviaua i 


failure rates used 




Mu 






„ ^ # v is a J 






1 UU 1 V 1U Ud 1 OCi-VXUC id LCu UDCU 




Dynamic 1 
Mean 


5.11 


4.73 


4.16 


3.44 


Variance 


0.94 


1.49 


1.91 


2.00 


S .Dev. 


0.97 


1.22 


1.38 


1.41 


S_ 

X 


0.04 


0.05 


0.06 


0.06 


% Up A/C 


0.85 


0.79 


0.69 


0.57 


Dynamic 2 
Mean 


4.94 


4 . 69 


4.07 


3.51 


Variance 


1.18 


1.38 


1.88 


1.99 


S . Dev . 


1.09 


1.18 


1.37 


1.43 


s_ 

>c 


0.05 


0.05 


0.06 


0.06 


% Up A/C 


0.82 


0.78 


0.68 


0.59 


Dynamic 3 
Mean 


5.64 


5.19 


4.40 


3.52 


Variance 


0.64 


1.33 


2.16 


1.95 


S . Dev. 


0.80 


1.15 


1.47 


1.39 


S_ 

X 


0.04 


0.05 


0.07 


0.06 


% Up A/C 


0.94 


0.87 


0.73 


0.59 


FIFO 

Mean 


5.21 


4.65 


3.89 


3.24 


Var iance 


1.18 


2.16 


2.86 


2.66 


S . Dev . 


1.09 


1.47 


1.69 


1.63 


S_ 

X 


0.05 


0.07 


0.08 


0.07 


% Up A/C 


0.87 


0.77 


0.65 


0.54 


FRP-Low 

Mean 


4.16 


3.85 


3.69 


3.25 


Variance 


4.67 


5.12 


4.28 


3.86 


S . Dev. 


2.16 


2.26 


2.07 


1.97 


S_ 

X 


0.09 


0.10 


0 . 09 


0.09 


% Up A/C 


0.69 


0.64 


0.62 


0.54 


FRP-H.iqJl 

Mean 


4.15 


3.88 


3.31 


3.18 


Variance 


3.42 


3.63 


3.53 


2.93 


S . Dev. 


1.85 


1.91 


1.53 


1.71 


S_ 

X 


0.08 


0.09 


0 . 09 


0.08 


% Up A/C 


0.69 


0.65 


0.55 


0.53 
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TABLE B-l 



MODEL STATISTICS: EXPECTED NUMBER OF UP 

AIRCRAFT VS SIMULATION CASE NUMBER 





Case 3a 


Case 3b 


Case 3c 


Case 


Planes 


12 


12 


12 


12 


K-Factor 


1 


0 


-1 


-2 


Parts/P lane 


9 


9 


9 


9 


Lambda 












inuiviuiiai failure 


idles useo 




Mu 




-Individual Service 


rates used 




Dynamic 1 
Mean 


11.89 


11 .12 


9 . 66 


7.42 


Var iance 


0.15 


1.09 


1.97 


3 . 52 


S . Dev . 


0 . 38 


1.05 


1 .40 


1.88 


S 


0.02 


0.05 


0.06 


0 .08 


% Up A/C 


0.99 


0.93 


0.80 


0.62 


Dynamic 2 
Mean 


11.87 


11.06 


9.69 


7.68 


Variance 


0 . 21 


1 .21 


1.99 


3.47 


S . De v . 


0 .46 


1.10 


1 .41 


1.86 


S 


0.02 


0.05 


0.06 


0.08 


% Up A/C 


0.99 


0.92 


0.81 


0 .64 


Dynamic 3 
Mean 


12.00 


11.96 


11.73 


9.84 


Var iance 


0.00 


0.14 


0.64 


4 . 43 


S . Dev . 


0.00 


0.38 


0.80 


2.10 


S_ 

X 


0.00 


0.02 


0.04 


0.09 


% Up A/C 


1.00 


0.99 


0.98 


0.82 


FIFO 

Mean 


11.99 


11.71 


10.49 


7.88 


Var iance 


0.01 


0 .43 


2.17 


5.99 


S . Dev . 


0.10 


0 . 66 


1.47 


2.45 


S_ 

x 


0.00 


0.03 


0.07 


0.11 


% Up A/C 


0.99 


0.98 


0.87 


0 . 66 


FRP-Low 

Mean 


10.97 


7.83 


4.04 


0.84 


Variance 


4.81 


12.27 


11.30 


2 .94 


S . Dev . 


2.19 


3.50 


3 . 36 


1.68 


S_ 

X 


0.09 


0 . 16 


0.15 


0.08 


% Up A/C 


0.91 


0 .65 


0 . 34 


0.07 


FRP-Hiah 

Mean 


11.12 


8.27 


4.53 


1.32 


Variance 


3.95 


9 . 89 


11 . 43 


4.16 


S . De v . 


1.99 


3.15 


3.38 


2.04 


S_ 

x 


0.09 


0.14 


0.15 


0.09 


% Up A/C 


0.93 


0.69 


0.38 


0.11 
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TABLE B-l 



MODEL STATISTICS: EXPECTED NUMBER OF UP 

AIRCRAFT VS SIMULATION CASE NUMBER 





Case 3e 


Planes 

K-Factor 

Parts/Plane 

Lambda 

Mu 


12 

-2.5 

9 




Dynamic 1 
Mean 

Variance 
S . Dev . 

S_ 

x 


6.15 
5 .03 
2.24 
0.07 


% Up A/C 


0.51 


Dynamic 1 

Mean 

Var iance 
S . Dev . 

S_ 

X 


6.56 

4.81 

2.19 

0.06 


% Up A/C 


0.54 


Dynamic 3 
Mean 

Var iance 
S . Dev . 

S_ 

X 


8.07 

6.74 

2.59 

0.12 


% Up A/C 


0.67 


FIFO 

Mean 

Var iance 
S . Dev . 

S_ 

X 


6.16 

8.55 

2.92 

0.08 


% Up A/C 


0.51 


FRP-Low 

Mean 

Var iance 
S . Dev . 

S_ 

X 


0.43 

1.43 

1.19 

0.05 


% Up A/C 


0.04 


FRP-Hiah 

Mean 

Var iance 
S . Dev . 

S_ 

x 


0.78 

2.13 

1.46 

0.05 


% Up A/C 


0.06 



Case 3f 

12 

-3.0 

9 

-Individual Failur 
Individual Service 



4 .46 

5 . 30 

2.30 
0.07 

0 .37 



5.23 

5.72 

2.39 

0.07 

0.44 



6.16 

6.49 

2.55 

0.11 

0.51 



8.13 

2.85 

0.08 

0.36 



0.38 

1.03 

1.01 

0.05 

0.03 



1.11 

0.05 

0.04 



Case 3a 


Case 3h 


12 


12 


-3.5 


-4 . 0 


9 


9 


rates used- 




ates used- 




2.75 


1.36 


4.58 


2 . 37 


2.14 


1.54 


0.09 


0.07 


0.23 


0.11 


3.31 


1.84 


5.24 


3.12 


2.29 


1.77 


0.10 


0.08 


0.28 


0.15 


4.20 


2.45 


5.82 


3.73 


2.41 


1.93 


0.11 


0.09 


0.35 


0.20 


2.53 


1.26 


5.01 


2.38 


2.24 


1.54 


0.10 


0.07 


0.21 


0.11 


0.38 


0.51 


1.46 


1.33 


1.21 


1.15 


0.05 


0.05 


0.03 


0.04 


0.51 


0.61 


1.25 


1.27 


1.12 


1.13 


0.05 


0.05 


0.04 


0.05 
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TABLE B-l 



MODEL STATISTICS: EXPECTED NUMBER OF UP 

AIRCRAFT VS SIMULATION CASE NUMBER 





Case 


4a 


Case 4b 


Case 4c 


Case 


4d Case 


Planes 


6 




6 


6 


6 


6 


Spares® 


12 




9 


6 


3 


1 


Parts/Plane 


9 




9 


9 


9 


9 


Lambda 




Failure rate 


all components = 


0.02 


Mu 




-Service rate 


all components = 


0 .84 


Dynamic 1 


Mean 


6.00 




6.00 


5.98 


5.28 


3.82 


Variance 


0 . 00 




0 . 00 


0.03 


1 . 09 


1.96 


S . Dev . 


0.00 




0.00 


0.18 


1 .04 


1.40 


S_ 


0.00 




0.00 


0.01 


0.05 


0.06 


% Up A/C 


1.00 




1.00 


0.99 


0.88 


0 .64 


Dynamic 2 


Mean 


6.00 




6.00 


5.97 


5.27 


3 .77 


Var iance 


0.00 




0.00 


0.05 


1.16 


2.13 


S .Dev. 


0.00 




0.00 


0.22 


1.08 


1.46 


S_ 

9€ 


0 . 00 




0.00 


0.01 


0.05 


0.07 


% Up A/C 


1.00 




1.00 


0.99 


0.88 


0.63 


Dynamic 3 


Mean 


6 . 00 




6.00 


5.96 


5.18 


3.63 


Variance 


0.00 




0.00 


0.06 


1.26 


2.19 


S . De v . 


0 . 00 




0 . 00 


0.24 


1.12 


1.48 


S_ 

9C 


0 . 00 




0 . 00 


0.01 


0.05 


0.07 


% Up A/C 


1.00 




1.00 


0.99 


0 . 86 


0.61 


FIFQ 

Mean 


6 . 00 




5.97 


5.78 


4.70 


3.33 


Variance 


0.00 




0 .06 


0.58 


2.32 


2.86 


S . Dev . 


0.00 




0.23 


0.76 


1.52 


1.69 


S_ 

yc 


0.00 




0.01 


0.03 


0.07 


0.08 


% Up A/C 


1.00 




0.99 


0.96 


0.78 


0.56 


FRP-Low 


Mean 


5.68 




5.09 


4 . 28 


3.65 


3.03 


Variance 


1.04 




2.78 


4.47 


4 . 49 


3.89 


S . Dev . 


1.02 




1.67 


2.11 


2.12 


1.97 


S_ 

yc 


0.05 




0 .07 


0.09 


0.09 


0 . 09 


% Up A/C 


0.95 




0.84 


0.71 


0.61 


0.51 


FRP-Hiah 


Mean 


5.68 




5.09 


4.28 


3.65 


3.03 


Var iance 


1.04 




2.78 


4.47 


4 . 49 


3.89 


S . De v . 


1.02 




1.67 


2.11 


2.12 


1.97 


S_ 

yc 


0.05 




0.07 


0.09 


0.09 


0.09 


% Up A/C 


0.95 




0.84 


0.71 


0 . 61 


0.51 



“Spares to Plane Ratio method used to calculate spares 
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TABLE B-l MODEL STATISTICS: EXPECTED NUMBER OF UP 

AIRCRAFT VS SIMULATION CASE NUMBER 

Case 5a Case 5b Case 5c Case 5d Case 5e 

Planes 12 12 12 12 12 

Spares 7 24 18 12 6 2 

Parts/Plane 99999 

Lambda Failure rate all components = 0.02 

Mu Service rate all components = 0.84 



Dynamic 1 
Mean 

Variance 
S . Dev. 

S_ 

SC 


12.00 

0.00 

0.00 

0.00 


11.92 

0.16 

0.39 

0.01 


% Up A/C 


1.00 


0.99 


Dynamic 2 
Mean 

Variance 
S . Dev . 

S 

sc 


12.00 

0.00 

0.00 

0.00 


11.94 

0.09 

0.29 

0.01 


% Up A/C 


1.00 


0.99 


Dynamic 3 
Mean 

Variance 
S . Dev . 

S_ 

sc 


12 . 00 
0.00 
0.00 
0.00 


11.91 

0.19 

0.44 

0.02 


% Up A/C 


1.00 


0.99 


FIFO 

Mean 

Variance 
S . Dev. 

S_ 

sc 


11.92 

0.25 

0.50 

0.02 


10.51 
4 . 57 
2 . 14 
0.10 


% Up A/C 


0.99 


0.88 



8.88 


3.15 


0.82 


3.86 


3.62 


0.89 


1.97 


1.90 


0.95 


0 . 09 


0.09 


0.04 


0.74 


0.26 


0.07 


8.95 


3.15 


0.81 


4.07 


3.75 


0.93 


2.02 


1.94 


0.96 


0.09 


0.09 


0.04 


0.75 


0.26 


0.07 


8.75 


2.94 


0.74 


4.01 


3.32 


0.72 


2.00 


1.82 


0.85 


0.09 


0.08 


0.04 


0.73 


0.25 


0.12 


5.61 


1.25 


0.66 


8.03 


2.07 


0.83 


2.83 


1.44 


0.91 


0.13 


0.06 


0.04 


0.47 


0.10 


0.06 



7 Spare to Plane Ratio method used for computing spares 
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TABLE B-l 



MODEL STATISTICS: EXPECTED NUMBER OF UP 

AIRCRAFT VS SIMULATION CASE NUMBER 





Case 6a 


Case 


6b 


Case 6c 


Case 6d 


Case 


Planes 


6 


6 




6 


6 


6 


Spares® 


12 


9 




6 


3 


1 


Parts/Plane 


9 


9 




9 


9 


9 


Lambda 


Individual 


component Failure rates 


used- 


Mu 


Individual 


component Service rates 


used- 


Dynamic 1 














Mean 


6.00 


6.00 




5.95 


5.31 


3.99 


Var iance 


0 . 00 


0.00 




0 . 09 


1.27 


2.05 


S . Dev. 


0.00 


0 . 00 




0.29 


1.13 


1.43 


S_ 

X 


0.00 


0.00 




0.01 


0.05 


0.06 


% Up A/C 


1.00 


1.00 




0.99 


0.89 


0 . 67 


Dynamic 2 














Mean 


5.99 


5.98 




5.81 


5.01 


3.78 


Var iance 


0.00 


0 . 07 




0 .54 


1.87 


2.64 


S . Dev. 


0.04 


0.26 




0.74 


1.37 


1.63 


S_ 

X 


0.00 


0.01 




0.03 


0.06 


0.07 


% Up A/C 


0 . 99 


0.99 




0.97 


0.84 


0.63 


Dynamic 3 














Mean 


6.00 


6.00 




5.95 


5.25 


3 . 84 


Var iance 


0.00 


0.00 




0.07 


1.36 


2 . 35 


S . Dev. 


0.00 


0.00 




0.26 


1.17 


1 . 53 


S_ 

X 


0.00 


0.00 




0.01 


0.05 


0.07 


% Up A/C 


1.00 


1.00 




0.99 


0.88 


0.64 


FIFO 

Mean 


5.97 


5.90 




5.50 


4.54 


3.37 


Variance 


0.07 


0.29 




1.36 


3.15 


3.49 


S . Dev . 


0.26 


0.55 




1.16 


1.78 


1.87 


S_ 

X 


0.01 


0.02 




0.05 


0.08 


0.08 


% Up A/C 


0 .99 


0 .98 




0.92 


0.76 


0.56 


FRP-Low 














Mean 


4.61 


4.28 




4.22 


4.04 


3.26 


Variance 


4.94 


5.45 




5.55 


5.09 


4.77 


S . Dev. 


2.22 


2.34 




2.36 


2.26 


2.18 


S_ 

x 


0.10 


0.10 




0.11 


0.10 


0.10 


% Up A/C 


0.77 


0.71 




0.70 


0.67 


0.27 


FRP-HiQh 














Mean 


5.63 


5.53 




5.10 


4 . 39 


3 . 50 


Var iance 


1 . 45 


1.72 




2.81 


3.83 


3.75 


S .Dev. 


1.21 


1.31 




1.68 


1.96 


1.94 


S_ 

X 


0.05 


0.06 




0.07 


0.06 


0.09 


% Up A/C 


0.94 


0.92 




0.85 


0.73 


0.58 



“Spare to Plane Ratio method used to calculate spares 
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TABLE B-2 RANDOM SEED EFFECTS: EXPECTED NUMBER OF UP 
AIRCRAFT AT MISSION TERMINATION* 





Case A 


Case B 


Case C 


Case D 


Planes 


12 


12 


12 


12 


Spares 10 


3 


3 


3 


3 


Parts/Plane 9 
Random 


9 


9 


9 


Seeds 


16989 


15989 


14989 


14704 




14997 


14999 


14731 


13331 


Dynamic 1 


13997 


13999 


16789 


15652 


Mean 


1.95 


1.87 


2.02 


1.87 


Variance 


2.87 


3.19 


3 . 40 


2 .72 


S . Dev. 


1.69 


1.78 


1.84 


1.65 


s . 11 

X 


0.08 


0.08 


0.08 


0.07 


95% Cl 12 


2.11,1.80 


2.04,1.72 


2.19,1.86 


2.02,1.73 


Dynamic 2 


Mean 


0.90 


0.83 


0.85 


0.95 


Var iance 


1.68 


1.80 


1 .76 


1.88 


S . Dev . 


1.29 


1.34 


1.32 


1.37 


S_ 

X 


0.06 


0.06 


0.06 


0.06 


95% Cl 


1.02,0.79 


0.95,0.72 


0.97,0.73 


1.07,0.83 


FIFO 


Mean 


0.76 


0.82 


0.75 


0.83 


Variance 


1.48 


1.52 


1.13 


1.40 


S . Dev . 


1.21 


1.23 


1.06 


1.18 


S_ 

X 


0.05 


0.06 


0.05 


0.05 


95% Cl 


0.87,0.66 


0.94,0.72 


0.85,0.66 


0.94,0.73 


FRP-Low 


Mean 


1.61 


1.82 


1.84 


1.76 


Var iance 


6.68 


7.26 


7.72 


7.60 


S . Dev . 


2.58 


2.69 


2.77 


2.75 


S_ 

X 


0.12 


0 . 12 


0.12 


0.12 


9 5% Cl 


1.84,1.38 


2.06,1.59 


2.09,1.60 


2.01,1.53 


FRP-Hiah 


Mean 


1.57 


1.65 


1.69 


1.68 


Var iance 


4.91 


4 . 55 


4.61 


4.79 


S . Dev . 


2.21 


2.13 


2.14 


2.18 


S_ 

X 


0.10 


0.10 


0.10 


0.09 


9 5% Cl 


1.76,1.38 


1.84,1.47 


1.88,1.51 


1.88,1.49 



*500 replications 

10 Equal numbers of spares for all type components. 
1X S- = { ( Var iance/500 ) } - * = Standard Error 

x 

12 95% Cl = Mean t S_ x 1.96 , table values are 
endpoints 
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TABLE 


B-3 


SPARE 


ALLOCATION 


OF COMPONENTS BY 










CASE 


NUMBER: K-STANDARD DEVIATION 


METHOD 




Case 


l 


Case lb 


Case lc 


Case 


Id 


K Factor 


+ 1.0 




0.0 


-0 . 5 


-1.0 




PART 1 


10 




7 


6 


5 




PART 2 


7 




5 


4 


3 




PART 3 


9 




6 


5 


4 




PART 4 


2 




1 


0 


0 




PART 5 


10 




7 


6 


4 




PART 6 


11 




8 


6 


4 




PART 7 


1 




0 


0 


0 




PART 8 


6 




4 


3 


2 




PART 9 


14 




11 


9 


8 






Cass. 


le 


Case If 


Case la 


Case 


lh 


K Factor 
PART 1 


-1. 

3 


5 


-2.0 

2 


-2 . 5 
0 


-3.0 

0 




PART 2 


2 




1 


0 


0 




PART 3 


3 




1 


0 


0 




PART 4 


0 




0 


0 


0 




PART 5 


3 




1 


0 


0 




PART 6 


4 




2 


1 


0 




PART 7 


0 




0 


0 


0 




PART 8 


1 




0 


0 


0 




PART 9 


6 




4 


3 


1 






Case 


2a 


Case 2b 


Case 2c 


Case 


2 <4 


K Factor 


+ 1.0 




0.0 


-1.0 


-1.5 




PART 1 


18 




14 


11 


9 




PART 2 


13 




10 


7 


5 




PART 3 


17 




13 


9 


8 




PART 4 


3 




2 


0 


0 




PART 5 


18 




14 


10 


8 




PART 6 


20 




16 


12 


10 




PART 7 


2 




1 


0 


0 




PART 8 


11 




8 


5 


4 




PART 9 


26 




22 


17 


15 






Casa. 


_2s. 


Case 2f 


Case 2a 


Case 


2h 


K Factor 


-2.0 




-2 . 5 


-3.0 


-3 . 5 




PART 1 


7 




5 


3 


1 




PART 2 


4 




2 


1 


0 




PART 3 


6 




4 


2 


0 




PART 4 


0 




0 


0 


0 




PART 5 


8 




5 


3 


1 




PART 6 


8 




6 


4 


2 




PART 7 


0 




0 


0 


0 




PART 8 


2 




1 


0 


0 




PART 9 


12 




10 


8 


5 
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TABLE B-3 MODEL PARAMETERS: SPARE ALLOCATION OF 

COMPONENTS VS SIMULATION CASE NUMBER 





Case 


3a 


Case 


3b 


Case 


3c 


Case 




K Factor 


+ 1.0 




0.0 




-1.0 




-2 . 0 




PARTI 


34 




29 




23 




18 




PART 2 


25 




21 




16 




12 




PART 3 


31 




26 




21 




16 




PART 4 


6 




4 




2 




0 




PART 5 


33 




28 




23 




17 




PART 6 


37 




31 




26 




20 




PART 7 


3 




1 




0 




0 




PART 8 


20 




16 




12 




8 




PART 9 


50 




43 




37 




30 






Case 


3e 


Case 


3f 


Case 


-3a 


Case 


3h 


kBFP 1 - 


-2 . 5 
15 




-3.0 

13 




-3.5 

10 




-4.0 

7 




PART 2 


9 




7 




5 




2 




PART 3 


13 




11 




8 




6 




PART 4 


0 




0 




0 




0 




PART 5 


15 




12 




9 




7 




PART 6 


17 




15 




12 




9 




PART 7 


0 




0 




0 




0 




PART 8 


6 




4 




2 




0 




PART 9 


27 




23 




20 




17 
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APPENDIX C 



TABLE C-l. FIFO ALGORITHM 
Key Variable Description; 

Planes: The total number of whole systems operational, in 
this case aircraft. 

TOTSPR: The initial number of spares per type part. 

TOTMSN: The total mission length in flight hours (120 Hrs.). 

TOTPRT: The total number of types of parts. 

TOTRUN: The total number of replications or runs (100). 

MSNTIM: The mission time running clock, ranges vary from 
0 to 120 hours (TOTMSN). 

LOWLVL: An Array which stores the lowest part stock level in 
the event of a given part failure. 

STOCK: An Array which stores the number of up and working 
components of a given part type (spares + planes). 

FAILT: A variable which holds the failure time generated by 
a uniform random variable and the inverse transform 
method for an exponential distribution. This is 
compared to the repair time to determine the 
scheduling of next event. 

COMPLT: A variable which holds the repair time, is generated 
like the failure time, and when compared to the 
latter determines the next scheduled event. 

SERVTM: A variable which is assigned a exponential value 
generated by the inverse transform method and is 
added to the current clock to determine the 
completion time of the next repair 

QPART: An array which stores component types as the wait in 
the queue. 

QSIZE: An array which stores current position a component 
holds in the queue. 

QTIME: An array which stores the time a part enters the 
queue . 

RHO: The traffic intensity equal to arrival rate divided by 
the service rate in tne queue 

SUMLAM: The sum of all the part failure rates. 

SRATE : the service rate ()i) in the queue. 

REPAIR: A counter to track when repairs occur. 

FLAG2 : A counter to track if stocks are depleted. 

RUN: A counter to track run number (0->100) 

Input: Failure rate parameters (Lambda(i)) for each 
component . 
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Repair rate parameters (SRATE(i)) for each component. 



Output : 



Step 1: 

Step 2: 
Step 3: 
Step 4 : 

Step 5: 
Step 6: 

Step 7: 
Step 9 

Step 10 



Mean number of up aircraft at missions end 

Variance and standard deviations of above means 

Distribution of up aircraft at missions termination 
(the number of runs with zero aircraft up, the 
number with two, and so forth). 

The average waiting time in the queue 

The average queue line 

Initialize simulation parameters: TOTRUN, TOTMSN, 
TOTPRT, SR ATE, etc. 



Read in repair rates SRATE(i) and failure rates 
LAMBDA( i ) 

Call RANDOM; Generate an array to store uniform 
random numbers, store in A( ),AA( ),AAA( ). 

Initialize a new mission run: MSNTIM, PLANES, 

SPARE ( I ) , COMPLT, FLAG2, 

STOCK = SPARES + PLANES 

SUMLAM = LAMBDA ( 1 ) + LAMBDA( 2 ) + ... + LAMBDA(N) 

RHO = SUMLAM * PLANES / SR ATE 13 

Verify at least one aircraft is up (FLAG2 = 1) 

Generate a failed part: 

FAILT = - In (U) / (SUMLAM X PLANES) 

This uses the inverse transform method with a 
uniform random number (U) from an IMSL routine 
LLRANDOM II to yield an exponential interarrival 
t ime . 

Scan stock levels (STOCK(i)) and verify at least one 
aircraft is up (FLAG2 = 1): 

If STOCK ( i ) < 1 (FLAG2 =0), then check to see if a 

repair has occurred which will create an up 
aircraft. If all aircraft are down, advance to next 
repair time or terminate run. 

Test to see if a failure or repair is the next 
event to occur: 

Is FAILT < COMPLT ? If so, MSNTIM = FAILT. go to 
step 11; else MSNTIM = COMPLT.go to step 10. Test 
both cases to see if mission has terminated: 

Is MSNTIM greater than TOTMSN ? If so, go to step 
18 . 

Call ALTQ; Enter gueue, complete repair and 
increment appropriate stock. Commence service on 



Not calculated when using multiple service rates. 
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the next part if components remain waiting 
service, go to step 7 . 



Step 11: Determine which part failed; calculate the 

individual probability a given part i failed: 

Probability Part i failed = LAMBDA( i ) /SUMLAM 

Then, generate a discrete CDF with the 

S robabilities. Pick a uniform random variable to 
etermine the type failure. Decrement the 
appropriate stock. 

Step 12: Call QUEUE: Enter queue with failed part. 

STEP 13: Check to see if server is available and the queue 

size is less than one: 

Is SERVR = 0 and QSIZE = 0 ? If so. enter service, 
go to step 15, else join queue, go to step 14. 

Step 14: Part joins queue at the end of the service 

line. Increment QSIZE, store QPART, QTIME, check 
QMAX versus QSIZE, go to step 7. 

Step 15: Enter service, store waiting time in queue 

(WAIT(i)), part type (QPART(i)), advance positions 
in the queue up one position. Calculate service 
t ime : 

SERVTM = - In (U) / SRATE(i) 

Calculate completion of service time (COMPLT): 
COMPLT = MSNTIM + SERVTM 



Step 16: Store the stock levels at the moment of a given 

failure (JFAIL) at a given time (MSNTIM). Store in 
LOWLVL the lowest stock level at that time. 



Step 17: If stocks are depleted, cease failures until a 

repair of the proper type occurs. FLAG2 = 0, Go to 
step 10. 

Step 18: Call ALTQ1 : Talley waiting time, queue size of 

those parts which remain in the queue at mission 
termination time MSNTIM 

Step 19: Is this run the last run: RUN i TOTRUN ? 

If not, go to step 5 

Step 20: Call STATS: compute statistics; 

Calculate mean up aircraft: MUPPLN = UPPLANE/RUN 

Calculate the variance and standard deviation of 
MUPPLN 



Step 21: Calculate distribution of up aircraft; talley the 

number of up aircraft less than one standard 
deviation from the mean, two standard deviations 
from the mean, or at the mean, etc (CNTR). Talley 
the number of runs with zero aircraft up, one 
aircraft up, and so forth (YCNTR). 
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Step 22: 

Step 23: 
Dynamic 1 

Step 14A: 

Step 14b: 

Step 14c: 

Dynamic 2 

Step 14a: 

Step 14b: 



Step 14c: 

Dynamic 3 
Step 14a: 

Step 14b: 



Step 14c: 

FRL-L(H) 



Calculate the average wait in line WAITT, average 
maximum qsize MQBAR, and average number of 
customers waiting in the queue. 

Print Statistics 



Model Modifications 14 : 



Count the number of parts of type i in the queue 
line. Are there equal numbers of parts ? If so, go 
to step 14c, else go to step 14b. 

Sort components by count, reorder position in line 
in favor of the most numerous component (to the 
front of the line). Return to step 7. 

Sort parts by arrival time in the line (FIFO). 
Return to step 7. 

Model Modifications: 

Count the number of parts of type i in the queue 
line ( COUNT ( i ) ) . 

Assign a priority to each type component i: 

RHO(i) = LAMBDA ( i ) / SRATE(i) 

LINEUP ( i ) = RHO(i) x COUNT(i) 

Sort components by LINEUP(i), reorder position in 
line in favor of the largest LINEUP(i). Return to 
step 7. 

Model Modifications: 



Count the number of operating parts of each 
component i (STOCK(i)). 

Assign a priority to each component type based on 
the count of current stock in step 14a: 

LINEUP ( i ) < LINEUP ( j ) iff 

STOCK ( i ) < STOCK ( j ) 



Sort the parts in line by LINEUP(i), 
iositions in favor of the smallest LI 
eturn to step 7. 



reorder 



NEUP ( i ) . 



Model Modifications: 



Step 2b: Call LINEUP; Read in the values of LAMBDA(i), 

assign a priority from one to n total type 
components based on the smallest (largest) 
LAMBDA(i). Return to step 3 

Step 14a: Sort queue line by priority values assigned in 
step 2b. Reorder position in line to favor the 
components with the smallest LAMBDA(i). Return to 
step 7. 



14 Modif ications apply to FIFO algorithm, steps follow 
sequentially from there 
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APPENDIX D 



Table D-l. FIFO Fortran Coding 



X INITIALIZATION OF VALUES 

INTEGER RUN, TOTRUN, TOTSPR, FPART , QMAXMNC 50 0 ) , 
1CNT1,CNT2,CNT3,CNT4, YCNT1 , YCNT2, YCNT3, YCNT4 , YCNT5 , YCNT6 , YCNT7 , 
1YCNT8,YCNT9,YCNT10,YCNT11,YCNT12,YCNT13 

REAL PLANES, L AMBDAC 10), SPAREC 10),STOCK(10),UPLAN( 500), LQ,K, 
1AC100000), LOWLVLC 500), TOTAL (10, 500), AAC 10 0000), CWAITC 1000 00), 
1PROBC 10) ,MSNTIM,MQBAR, QTC 50 0 , 120), FACTORC 10) , 

1PPLANE, AAA (60000), MU PLN, SR AT EC 10) 

COMMON/QUEU/STOCK,MSNTIM,L RUN, QMAXMN, LI NEK, CWAIT, CUSTMR, REPARE, 
1SYSARR, TZZ, QT, SRATE, TOTMSN, COMPLT, /STAT/MUPLN, VRPLN, SDPLN, CNT1 , 
1CNT2,CNT3,CNT4, YCNT1, YCNT2, YCNT3, YCNT4 , YCNT5, YCNT6 , YCNT7 , 
1YCNT3, YCNT9,YCNT10,YCNT11,YCNT12,YCNT13, 

1MQBAR, WAITT, LQ 

DATA JFAIL/0/,1/0/, II/0/,RUN/1/,TOTRUN/500/, 

1TOTPRT/9/, FLAG2/l/,K/-3.5/ 

TOTMSN = 120 
TOTSPR = 01 
PLANES =06.0 

x READ IN FAILURE AND REPAIR RATES 

DO 1 J = 1 , TOTPRT 

READ (2,2) LAMBDA(J) 

2 FORMAT ( F8 . 5 ) 

1 CONTINUE 

DO 3 J = 1, TOTPRT 

READ (3,8) SRATE(J) 

8 FORMAT ( F8 . 5) 

3 CONTINUE 

SUMLAM = 0 

DO 7 J = 1, TOTPRT 

SUMLAM = SUMLAM + LAMBDA(J) 

7 CONTINUE 

DO 9 I = 1, TOTPRT 

FACTORC I ) = PLANES * LAMBDA(I) * TOTMSN 

SPARE ( I ) = NINT ( FACTORC I ) + K * SQRTC FACTORC I )) ) 

IFC SPAREC I ) .LT. 0) SPARECI) = 0 

9 CONTINUE 

x CALL LLRAND II, GENERATE RANDOM NUMBER ARRAYS 

CALL RANDOMC A, AA, AAA) 
x INITIALIZE A NEW MISSION 

5 MSNTIM = 0.0 

COMPLT = 120 

J FAIL = 0 
PPLANE = PLANES 

XXX = SUMLAM * PLANES 

L RUN = RUN 

FLAG2 = 1 

FLAG3 = 1 

DO 6 J = 1, TOTPRT 

STOCKCJ) = SPAREC J ) + PLANES 

6 CONTINUE 
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REPARE = 0 



10 CONTINUE 

IFCFLAG2 .GE. 1) THEN 

* GENERATE ARRIVALS OF FAILED PARTS 

1 = 1 + 1 

Z = I 
U = A ( I ) 

T = -LOG( U ) / C SUMLAM * PLANES) 

FAILT = T + MSNTIM 
ENDIF 

^ TEST TO SEE IF A FAILURE HAS OCCURED 

20 CONTINUE 

I F( REPARE .GE. 1 .AND. FLAG2 .LT.l) THEN 
REPARE = 0 
DO 25 J = 1 , TOTPRT 

I F( STOCK ( J ) .LT. 1) THEN 
FLAG2 = 0 
FAILT = TOTMSN 
GO TO 23 

ELSE 

FLAG2 = 1 

ENDIF 

25 CONTINUE 
GO TO 10 
ENDIF 

REPARE = 0 

DO 26 J = 1, TOTPRT 

IF(STOCKCJ) .LT. 1) THEN 
FLAG2 = 0 
FAILT = TOTMSN 
GO TO 23 

ELSE 

FLAG2 = 1 

ENDIF 

26 CONTINUE 

* INCREMENT TIME STEP 

X VERIFY MISSION LENGTH 

23 IFCFAILT .LT. COMPLT .AND. FLAG2 .GE. 1)THEN 

MSNTIM = FAILT 

IF( MSNTIM .GE. TOTMSN) GO TO 22 
GO TO 30 
ELSE 

MSNTIM = COMPLT 

IFCMSNTIM .GE. TOTMSN) GO TO 22 
GO TO 21 
ENDIF 

21 CALL ALTQ 

29 GO TO 20 

x CALCULATE WHICH COMPONENT FAILED 

30 CONTINUE 

DO 31 J = 1 , TOTPRT 

PROBCJ) = LAMBDA ( J ) / SUMLAM 

31 CONTINUE 

X CALL A RANDOM NUMBER, USE LAMBDA’S TO BUILD CDF, 

X AND DETERMINE WHICH PART FAILED 



11 = 11+1 
ZZ = II 



32 

33 

X 

34 

x 

40 

41 

42 

x 

51 

x 

22 

70 

72 

71 

x 



X = AA( II ) 

XL NR = 0.0 

DO 32 J = h TOTPRT 

XUPPR = XL HR + PROB(J) 

I F( X .GT. XL HR .AND. X .LE. 
STOCK(J) = STOCK(J) - 
FPART = J 
GO TO 33 



XUPPR) THEN 

1 



ELSE 

XL HR = XUPPR 



ENDIF 

CONTINUE 



CALL QUEUEC FPART, AAA, *34, *29 ,*70) 



CHECK STOCK LEVELS AND BRANCH 



CONTINUE 

ARE STOCKS LEVELS = 0? 

CONTINUE 

DO 41 J = 1 , TOTPRT 

IFCSTOCK(J) .GE. 1) THEN 
FLAG2 = 1 

ELSE 

FLAG2 = 0 
FAILT = TOTMSN 
GO TO 42 

ENDIF 

CONTINUE 

JFAIL = JFAIL + 1 

RECORD STOCK LEVELS FOR A GIVEN FAILURE 

DO 51 J = 1, TOTPRT 

TOTALU, JFAIL) = STOCK(J) 

CONTINUE 

GO TO 10 



MISSION TERMINATION 
CALL ALTQ1 
CONTINUE 

DO 71 JJ = 1, JFAIL 
L OWL VL ( J J ) = 15 
DO 72 LL = 1, TOTPRT 

IFCLOWLVLCJJ).GT.TOTALCLL,JJ)) LOWLVL(JJ) = TOTALCLL, JJ) 

CONTINUE 

CONTINUE 

IFCLONLVLC JFAIL) .LE. PLANES ) THEN 
UPLANC RUN) = LOWLVLC JFAIL) 

ELSE 

UPLANC RUN) = PLANES 
ENDIF 



PLOT AND DISPLAY OUTPUT 
I F( RUN .GE. TOTRUN)THEN 

HR I TEC 06 , X) 'CONTINUOUS TIME, CONSTANT FAILURE RATE, FIFO 1 
HRITEC 06 , X) ‘MISSION TIME = ' , MSNTIM, 'JFAIL 1 , JFAIL 
HRITEC06,*) 'INITIAL PLANES = ' , PL ANES , ' I NITI AL SPARE FACT OR 1 , K 
HRITEC 06 , x ) 'RUN NUMBER * , RUN, * CUMULATIVE FAIURES CI)= f ,ZZ 



S3 



tt tt 



* 

* 

200 

* 

102 

103 



104 

199 



Z, ZZ , TZZ TRACK ARRAY INPUTS TO VERIFY THEIR DIMENSIONS ARE 

NOT EXCEEDED 

WRITEC06,*) 'TOT PARTS 3 ' , TOTPRT , 'Z 3 ',Z, 'ZZZ',TZZ 
WRITEC06,*) 1 L AMBDA= ' , XXX, • SERVICE RATE 3 VARIABLE' 
WRITEC06,*) 'TRAFFIC INTENSITY RHO= VARIABLE' 

ENDIF 

TEST RUN NUMBER TO STOP SIMULATION 

I FC RUN .GE. TOTRUN) GO TO 200 
RUN = RUN + 1 
GO TO 5 

COMPUTE SAMPLE STATISTICS 



CONTINUE 

CALL STATS ( CUSTMR, L RUN, UPL AN, QMAXMN, LINEK, CWAIT, 
1TOTMSN, SYSARR, QT) 

WRITE OUTPUT OF STATISTICS 

WRITE( 06,102) MUPLN, VRPLN, SDPLN 

FORMATC1X, 'EXPECTED NUMBER OF UP PLANES =»,F11.8, 
1/, IX, 'VARIANCE OF UP PLANES =',F15.4, 



1/, IX, 'STD. DEV. UP PLANES =',F15.4) 

WRI TEC 06,103) CNT1, 

1CNT2,CNT3,CNT4,RUN,T0TMSN, 

1,YCNT1,YCNT2,YCNT3,YCNT4,YCNT5,YCNT6,YCNT7, 
1YCNT8,YCNT9,YCNT10,YCNT11,YCNT12,YCNT13 
FORMATC1X, 'NUMBER OF RUNS UP PLANES <= -1 SIGMA 



,14, 



i/,ix, 


'NUMBER 


OF 


RUNS 


UP 


PLANES <= 


MEAN =',IA, 


1/,1X, 


•NUMBER 


OF 


RUNS 


UP 


PLANES <= 


+1 SIGMA =*,I4 


1/,1X, 


•NUMBER 


OF 


RUNS 


UP 


PLANES <= 


+2 SIGMA =’,I4 


1/,1X, 


'TOTAL RUNS =' , 


15, 


t 


MISSION 


TIME=',F5.0, 


1/,1X, 


'NUMBER 


OF 


RUNS 


W/ 


0 


PLANES 


UP = 1 , I 5, 


l/,ix. 


'NUMBER 


OF 


RUNS 


w/ 


1 


PLANE 


UP =',I5, 


1/,1X, 


•NUMBER 


OF 


RUNS 


W/ 


2 


PLANE 


UP , 15, 


1/,1X, 


'NUMBER 


OF 


RUNS 


w 


3 


PLANE 


UP = 1 , I 5 , 


l/.ix. 


'NUMBER 


OF 


RUNS 


w 


4 


PLANE 


UP = ’ , 15, 


1/,1X, 


'NUMBER 


OF 


RUNS 


w/ 


5 


PLANE 


UP = ’ , 15, 


1/, IX, 


'NUMBER 


OF 


RUNS 


w/ 


6 


PLANE 


UP = ’ , 15, 


1/,1X, 


•NUMBER 


OF 


RUNS 


w 


7 


PLANE 


UP -',15, 


1/,1X, 


'NUMBER 


OF 


RUNS 


w 


8 


PLANE 


UP = ' , 15, 


1/,1X, 


'NUMBER 


OF 


RUNS 


w 


9 


PLANE 


UP = ' , I 5, 


1/,1X, 


•NUMBER 


OF 


RUNS 


W/10 


PLANE 


UP = ' , I 5, 


1/,1X, 


'NUMBER 


OF 


RUNS 


W/ll 


PLANE 


UP = ' , 15, 


l/,ix. 


'NUMBER 


OF 


RUNS 


W/ 12 


PLANE 


UP = ' , 15) 


WRITEC 06,104) 


MQBAR,I 


WAITT, 


LQ 





FORMAT (IX, ' AVE MAX QUEUE SIZE =',F5.0, 

1/ , IX, ' AVE WAIT IN THE QUEUE =',F10.4, 

1/, IX, 'AVE NUMBER OF CUSTOMERS WAITING IN Q=',F10.4) 



DO 199 I = 1, TOTPRT 

PRINT * , ' SPARES TYPEC ' , I , ' ) = ' , SPAREC I ) 

CONTINUE 

PRINT X, ' K= * , K 

STOP 

END 






************* X ************************************** XX**************** 



SUBROUTINE QUEUEC FPT , QA,* , * , * ) 

SAVE QSIZE, QPART, LINEUP, SERVR, SERVPT,SERVTM,COMPLT, QMAX, RUN, II, 

1 ARRIVE, WAIT, QTIME 

INTEGER QSIZE, SPART , QPART ( 200 ) , SERVPT , 

1QMAXC 500), QRUN, RUN,CNTR, FPT 

REAL QA( 6 0000 ) , QSTOCKC 10 ) , QTC 500,1 20 ),WAIT(1 0000 0),QTIME< 200 ) , 
lSRATE(lO) 

COMMON /QUEU/QSTOCK,QMSNTM, QRUN, QMAX, CNTR, WAIT, QCUST, REPAIR, 
1SYSARR, ZZZ, QT, SRATE, TOTMSN, COMPLT 

DATA SERVR/0/,QSIZE/0/, SERVPT/O/, I J/0/ , RUN/O/ , I I/O/, 

1 FLAG/ 0/ 

************************************************************************ 

* PARTS ARRIVE 

I F ( QRUN .LT. 2 .AND. FLAG .LT. 1) THEN 
FLAG = 1 
QCUST = 0 
SYSARR = 0 
DEPART = 0 
COMPLT = TOTMSN 
ENDIF 

CLOCK = QMSNTM 
ARRIVE = 1 
SYSARR = SYSARR + 1 

* INITIALIZE QUEUE FOR A NEW RUN 
I F( QRUN .GT. RUN) THEN 

QMAX(QRUN) = 0 
QSIZE = 0 
RUN = QRUN 
SERVR = 0 
COMPLT = TOTMSN 
ENDIF 



* IS THE SERVER BUSY? 

I F( SERVR .LT. 1 .AND. QSIZE .LT. 1) THEN 
QPART ( 1 ) = FPT 
GO TO 10 
ENDIF 
GO TO 30 



* CHECK SERVER AVAI LABL E, QUEUE, AND SERVICE TIME 

ENTRY ALTQ 
ARRIVE = 0 
CLOCK = QMSNTM 

* INITIALIZE QUEUE FOR A NEW RUN 
I F( QRUN .GT. RUN) THEN 

QMAX(QRUN) = 0 
QSIZE = 0 
RUN = QRUN 
SERVR = 0 
COMPLT = TOTMSN 
ENDIF 

1 I F ( SERVR .LT. 1 .AND. QSIZE .GE. 1) GO TO 10 

GO TO 100 
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* * 



X 



THE QUEUES 



30 QCUST = QCUST + 1 

IFCQSIZE .LT. 1) GO TO 40 
QSIZE = QSIZE + 1 
QPART ( QSIZE) = FPT 
QTIMEC QSIZE ) = CLOCK 

IFCQMAXCQRUN) .LT. QSIZE) QMAXCQRUN) = QSIZE 
GO TO 1 

40 CONTINUE 

QSIZE = 1 
QPART ( 1 ) = FPT 
QTIMEC QSIZE ) = CLOCK 

I F ( QMAXC QRUN ) .LT. QSIZE) QMAXC QRUN) = 1 
GO TO 1 

X THE SERVICE AREA 

10 CONTINUE 

SPART = QPART ( 1 ) 

X STORE WAITING TIMES FOR QUEUE 

IFCQSIZE .GT. 0) THEN 
II = II + 1 

WA ITCH) = CLOCK - QTIMEC 1 ) 

X STORE NUMBER OF PEOPLE IN LINE IN VARIABLE CNTR 

ENDIF 

CNTR = II 

IFCQSIZE .GT. 1) THEN 
PAIRS = QSIZE - 1 

DO 50 I = 1, PAIRS 

QPART ( I ) = QPART C 1 + 1 ) 

QTIMEC I ) = QTIMEC 1+1 ) 

50 CONTINUE 

QPART C QSIZE) = 0 
QTIMEC QSIZE) = 0 
ENDIF 

IFCQSIZE .GT. 0) QSIZE = QSIZE -1 



SERVR = 1 

SERVPT = SERVPT + 1 

x SERVICE IS EXPONENTIAL 

IJ = IJ + 1 
XX = QACIJ) 

ZZZ = IJ 

SERVTM = -LOGCXX) / SRATEC SPART ) 

COMPLT = SERVTM + CLOCK 

CHECK TO SEE IF SERVICE HAS COMPLETED 
IF SO INCREMENT PART 

100 IFC COMPLT .LE. CLOCK) THEN 
COMPLT = TOTMSN 
REPAIR = 1 
SERVR = 0 

DEPART = DEPART + 1 
QSTOCKC SPART ) = QSTOCK ( SPART ) + 1 
ENDIF 



-,b 



* RETURN TO APPROPRIATE SEGMENT OF MAIN PROGRAM 

101 IFCREPAIR ,EQ. 1 .AND. QSIZE . GT . 0 ) THEN 
REPAIR = 2 
GO TO 10 
ENDIF 

QT (RUN, CLOCK ) = QSIZE 
I F ( ARRIVE .GE. 1) RETURN 1 
RETURN 2 

x COUNT RESIDUAL PARTS IN QUEUE AT MISSION’S END 

ENTRY ALTQ1 

I F( QSIZE .GT. 0) THEN 

DO 129 I = 1 , QSIZE 
11 = 11+1 
CNTR = II 

WAITCII) = TOTMSN - QTIME(I) 

129 CONTINUE 

ENDIF 
RETURN 3 



END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE STATS C QCUST , RRUN, UPPLAN, QMAX, LINE, SWAIT , 

1 TOTMSN, SYSARR, QT ) 

COMMON/ ST AT/MUPP LN, VARPLN, S DP LAN, CNTR 1 , CNTR2, 

1CNTR3, CNTR4, YCNTR1 , YCNTR2, YCNTR3, YCNTR4, YCNTR5, YCNTR6, YCNTR7 , 
1YCNTR8,YCNTR9,YCNT10,YCNT11,YCNT12,YCNT13, 

1MXQBAR, WAIT , LQ 

INTEGER RRUN , QMAXC 500 ) , 

1CNTR1 , CNTR2,CNTR3,CNTR4, YCNTR1, YCNTR2, YCNTR3, YCNTR4, 

1YCNTR5 , YCNTR6 , YCNTR7 , 

1YCNTR8 , YCNTR9, YCNT10, YCNT11 , YCNT12, YCNT13 

REAL MUPPLN, PSUMC 500 ),UPPLAN( 500 ),MXQBAR, SWAIT Cl 00 000), 
1T0TMSN,LQ,QT(500,120) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

WAIT = 0 
CNTR1 = 0 
CNTR2 = 0 
CNTR3 = 0 
CNTR4 = 0 
YCNTR1 = 0 
YCNTR2 = 0 
YCNTR3 = 0 
YCNTRA = 0 
YCNTR5 = 0 
YCNTR6 = 0 
YCNTR7 = 0 
YCNTRS = 0 
YCNTR9 = 0 
YCNT10= 0 
YCNT 1 1 = 0 
YCNT12= 0 
YCNT13= 0 
MUPPLN = 0 
VARPLN = 0 
SDPLAN = 0 
MX Q BAR = 0 
LQ = 0 
UPLANE = 0 
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* >K 



SSPIAN = 0 
Q = 0 



x MEAN CALCULATIONS 

SUM UP A/C AT MSN TERMINATION FOR ALL RUNS 
SUM UP QUEUE MAX * S FOR EACH RUN 

DO 104 J = 1 , RRUN 

UPLANE = UPLANE + UPPLANCJ) 

Q = Q + QMAX(J) 

104 CONTINUE 

IF( UPLANE .GT. 0) THEN 
MUPPLN = UPLANE/RRUN 
ELSE 

MUPPLN = 0 
ENDIF 

DO 105 J = 1 > RRUN 

PSUM(J) = (MUPPLN - UPPLAN(J))x*2 
SSPLAN = SSPLAN + PSUMCJ) 

105 CONTINUE 

x MEAN MAX Q SIZE FOR ALL RUNS 

MXQBAR = Q / RRUN 



PRINT x, 1 LINE 1 , LINE, 1 MXQBAR f , MXQBAR , *Q f ,Q 

x CALCULATE THE AVE. WAIT IN QUEUE 

DO 108 IL = INLINE 

WAIT = WAIT + SWA IT ( IL ) 

108 CONTINUE 

WAIT = WAIT / SYSARR 

x ASSOCIATED VARIANCE/STD DEV CALC'S 

I F( SSPLAN .GT. 0 .AND. RRUN .GT.l) THEN 
VARPLN = SSPLAN / (RRUN - 1) 

SDPLAN = VARPLNX* . 5 
ELSE 

VARPLN = 0.0 
SDPLAN = 0.0 
ENDIF 

XI = MUPPLN - SDPLAN 
X2 = MUPPLN 
X3 = MUPPLN + SDPLAN 
X4 = MUPPLN + (2XSDPLAN) 

DO 106 J = 1 , RRUN 

IF( UPPLAN( J ) .LE. XI) CNTR1 = CNTR1 + 1 
IF(UPPLAN(J) . LE.X2 .AND. UPPLAN( J ) . GT . XI ) CNTR2 
IF(UPPLAN( J) .LE.X3 .AND. UPPLAN( J ) . GT . X2) CNTR3 
IF(UPPLAN( J) .LE.X4 .AND. UPPL AN( J ) . GT . X3 ) CNTR4 



106 CONTINUE 



CNTR2 + 1 
CNTR3 + 1 
CNTR4 + 1 
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107 



DO 107 J = 1 , RRUN 



IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
IFCUPPLANC J 
CONTINUE 



) . LE. 0) YCN 
) . L E . 1 .AND. 

) . L E . 2 .AND. 

) . L E . 3 .AND. 

) . L E . A .AND. 

) . L E . 5 .AND. 

) . L E . 6 .AND. 

) . L E . 7 .AND. 

) . L E . 8 .AND. 

) . L E . 9 .AND. 

) .LE.10 .AND. 

) . LE.ll .AND. 

) . L E . 12 . AND . 



TR1 = YCNTR1 
UPPLANC J) .GT. 
UPPLANC J ) .GT , 
UPPLANC J) .GT, 
UPPLANC J ) .GT . 
UPPLANC J ) .GT. 
UPPLANC J ) .GT. 
UPPLANC J) .GT. 
UPPLANC J ) .GT. 
UPPLANC J) .GT. 
UPPLANC J) .GT. 
UPPLANC J ) .GT. 
UPPLANC J) .GT. 



+ 1 

0) YCNTR2 = 

1) YCN TR3 = 

2) YCNTR4 = 

3) YCNTR5 = 

4) YCNTR6 = 

5) YCNTR7 = 

6) YCNTR3 = 

7) YCNTR9 = 

8) Y C N T 1 0 = 

9) YCNT 1 1 = 

10) YCNT12= 

11) YCNT1 3 : 



YCNTR2 
YCNTR3 
YCNTR4 
YCNTR5 
YCNTR6 
YCNTR7 
YCNTR8 
YCNTR9 
YCNT10+ 1 
YCNT11+ 1 
YCNT12+ 1 
YCNT13+ 1 



+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 



1 

1 

1 

1 

1 

1 

1 

1 



TEMP = 0 

DO 208 I = 1 , RRUN 

DO 209 J =1 f TOTMSN 

TEMP = TEMP + QTC I , J ) 
209 CONTINUE 

208 CONTINUE 

LQ = TEMP / C RRUN x TOTMSN) 

RETURN 

END 



XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX' 
SUBROUTINE RANDOMCRA, RAA,RAAA) 

REAL RAC1Q0000), RAAC 100 000 ), RAAAC 6 0000 ) 

N = 100000 

NN = 60000 

IX = 15989 

IIX = 14999 

IIIX = 13999 

CALL LRNDCIX,RA,N,1, 0) 

CALL LRNDCIIX,RAA,N,1,0) 

CALL LRNDCIIIX,RAAA,NN,1,0) 

RETURN 

END 
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A. 



Dynamic Modifications 
THE QUEUES 



x 



30 QCUST = QCUST + 1 

I F( QSIZE AT. 1) GO TO 40 
X STORE DATA OH NEWLY ARRIVED PART 

QSIZE = QSIZE + 1 
QPART (QSIZE) = FPT 
QTIME(QSIZE) = CLOCK 

IF(QMAX(QRUN) .LT. QSIZE) QMAX(QRUN) = QSIZE 

* ASSIGN Q PRIORITY BASED ON NUMBERS IN THE Q 

X INITIALIZE AND COUNT NUMBERS OF PARTS IN Q 

DO 31 I = 1, TOTPRT 
QCNTR(I) = 0 

31 CONTINUE 

DO 32 I = 1, TOTPRT 
DO 321 J = 1, QSIZE 

I F( I .EQ. QPART ( J ) ) THEN 

QCNTR(I) = QCNTR(I) + 1 

ENDIF 

321 CONTINUE 

32 CONTINUE 

DO 341 I = 1 , TOTPRT-1 

IFCQCNTR(I) .EQ. QCNTR( 1+1 ) ) THEN 
FFLAG = 0 

ELSE 

FFLAG = 1 
GO TO 37 

ENDIF 
341 CONTINUE 



x 

37 



332 

331 

x 

x 



SORT AND ASSIGN PRIORTIES OF SERVICE 
IF( FFLAG .GE. 1) THEN 



DO 331 I = 1, TOTPRT 

DO 332 J = 1, QSIZE 

I F( QPART ( J ) .EQ. I) THEN 
LINEUP(J) = QCNTRCI) 



ENDIF 

CONTINUE 

CONTINUE 



ASSIGN PRIORITY TO THE LINEUP WAITING IN THE Q 

SORT AND ASSIGN ACTUAL POSITIONS IN Q BASED ON PRIORITY 



PAIRS = QSIZE - 1 
DONE = 1 

35 IF( DONE .EQ. 1) THEN 

DONE = 0 

DO 351 I = 1 , PAIRS 

IFCLINEUP(I) .LT. LINEUP(I+1)) THEN 
TEMPL1 = LINEUP(I) 

TEMPL2 = QPART ( I ) 

TEMPL3 = QTIME(I) 

LINEUP! I ) = LINEUP(I+1) 

QPART ( I ) = QPART ( 1 + 1 ) 

QTIMECI) = QTIMEU + l) 

LINEUP! 1+1) = TEMPL1 
QPART ! 1+1 ) = TEMPL2 
QTIME! 1+1 ) = TEMPL3 
DONE = 1 

ENDIF 



60 



351 



36 



361 



40 

B. 

x 

30 

X 

* 

x 

31 



321 

32 



CONTINUE 

PAIRS - PAIRS - 1 
GO TO 35 

ENDIF 
END IF 



IFCFFLAG .LE. 0) THEN 
PPAIRS = QSIZE - 1 
PDONE = 1 

IFCPDONE .EQ. 1) THEN 
PDONE = 0 

DO 361 I = 1, PPAIRS 

IFCQTIMECI) .GT. QTIMECI+1)) 
PTEMP2 = QPARTCI) 

PTEMP3 = QTIMECI) 
QPARTCI) = QPARTCI+l) 
QTIMECI) = QTIMECI+1) 
QPARTCI+1) = PTEMP2 
QTIMECI+1) = PTEMP3 
PDONE = 1 



THEN 



ENDIF 

CONTINUE 

PPAIRS = PPAIRS - 1 
GO TO 36 



ENDIF 

ENDIF 



GO TO 1 

CONTINUE 
QSIZE = 1 
QPARTCI) = FPT 
QTIMEC QSIZE) = CLOCK 

IFC QMAXC QRUN) .IT. QSIZE) QMAXC QRUN) = 1 
GO TO 1 



Dynamic 2 Modifications 



THE QUEUES 

QCUST = QCUST + 1 

IFC QSIZE .LT. 1) GO TO 40 

STORE DATA ON NEWLY ARRIVED PART 
QSIZE = QSIZE + 1 
QPART C QSIZE) = FPT 
QTIMEC QSIZE) = CLOCK 

I FC QMAXC QRUN) .LT. QSIZE) QMAXC QRUN) = QSIZE 

ASSIGN Q PRIORITY BASED ON NUMBERS IN THE Q 

INITIALIZE AND COUNT NUMBERS OF PARTS IN Q 
DO 31 I = 1 , TOTPRT 
QCNTRCI) = 0 
CONTINUE 

DO 32 I = 1 , TOTPRT 
DO 321 J = 1 , QSIZE 

I FC I .EQ. QPART C J ) ) THEN 

QCNTRCI) = QCNTRCI) + 1 

ENDIF 

CONTINUE 

CONTINUE 
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SORT AND ASSIGN PRIORTIES OF SERVICE 



x 

DO 331 I = 1 , TOTPRT 

DO 332 J = 1 , QSIZE 

I F( QPART ( J ) ,EQ. I) THEN 

LINEUP(J) = QCNTRCI)xRHOCI) 

ENDIF 

332 CONTINUE 

331 CONTINUE 

X ASSIGN PRIORITY TO THE LINEUP WAITING IN THE Q 

x SORT AND ASSIGN ACTUAL POSITIONS IN Q BASED ON PRIORITY 

PAIRS = QSIZE - 1 

DONE = 1 

35 IF( DONE .EQ. 1) THEN 

DONE = 0 

DO 351 I = 1, PAIRS 

IF( LINEUPC I) .LT. LINEUPCI+1)) THEN 
TEMPL1 = LINEUP(I) 

TEHPL2 = QPART ( I ) 

TEMPL3 = QTIME(I) 

LINEUPC I ) = LINEUPC 1+1 ) 

QPART Cl) = QPART C 1 + 1 ) 

QTIMEC I ) = QTIMECI+1) 

LINEUPCI+1 ) = TEMPL1 
QPART C 1 + 1 ) = TEMPL2 
QTIMECI+1) = TEMPL3 
DONE = 1 

ENDIF 

351 CONTINUE 

PAIRS = PAIRS - 1 
GO TO 35 

ENDIF 
GO TO 1 

40 CONTINUE 
QSIZE = 1 
QPART Cl) = FPT 
QTIMEC QSIZE) = CLOCK 

IFCQMAXCQRUN) .LT. QSIZE) QMAXCQRUN)= 1 
GO TO 1 



C. Dynamic 3 Modifications 



X THE QUEUES 

30 QCUST = QCUST + 1 

IFC QSIZE .LT. 1) GO TO 40 
X STORE DATA ON NEWLY ARRIVED PART 

QSIZE = QSIZE + 1 
QPART C QSIZE) = FPT 
QTIMECQSIZE) = CLOCK 

IFCQMAXCQRUN) .LT. QSIZE) QMAXCQRUN) = QSIZE 
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COUNT NUMBERS OF OPERATING PARTS IN STOCK 



DO 341 I = 1 , T OTFRT- 1 

IF(QSTOCKCI) .EQ. QSTOCKCI+1)) THEN 
FFLAG = 0 

ELSE 

FFLAG = 1 
GO TO 37 

ENDIF 
341 CONTINUE 

^ SORT AND ASSIGN PRIORTIES OF SERVICE 

37 I F C FFLAG .GE. 1) THEN 

DO 331 I = 1 , TOTPRT 

DO 332 J = 1 , QSIZE 

IFCQPART(J) .EQ. I) THEN 

LINEUP(J) = QSTOCK(I) 

ENDIF 

332 CONTINUE 

331 CONTINUE 



x 

X 

35 



ASSIGN PRIORITY TO THE LINEUP WAITING IN THE Q 

SORT AND ASSIGN ACTUAL POSITIONS IN Q BASED ON PRIORITY 

PAIRS = QSIZE - 1 

DONE = 1 

I F ( DONE .EQ. 1) THEN 
DONE = 0 

DO 351 I = 1, PAIRS 

I F( LINEUPC I ) .GT. LINEUPCI+1)) THEN 
TEMPL1 = LINEUPCI) 

TEMPL2 = QPARTCI) 

TEMPL3 = QTIMECI) 

LINEUPCI) = LINEUPC 1+1 ) 

QPARTCI) = QPART C 1 + 1 ) 

QTIMECI) = QTIMEC 1+1 ) 

LINEUPCI+l) = TEMPL1 



351 



ENDIF 
ENDIF 



QPARTCI+1) 
QTIMEC 1+1 ) 
DONE = 1 

ENDIF 

CONTINUE 

PAIRS = PAIRS - 1 
GO TO 35 



TEMPL2 

TEMPL3 



IFC FFLAG .LE. 0) THEN 
PPAIRS = QSIZE - 1 
PDONE = 1 

36 IFC PDONE . EQ . 1) THEN 

PDONE = 0 

DO 361 I = 1 , PPAIRS 

IFC QTIMEC I ) .GT. QTIMECI + 1)) 
PTEMP2 = QPARTCI) 

PTEMP3 = QTIMECI) 
QPARTCI) = QPARTCI+1) 
QTIMECI) = QTIMECI+1) 
QPARTCI+1) = PTEMP2 
QTIMECI+1) = PTEMP3 
PDONE = 1 

ENDIF 

361 CONTINUE 



THEN 
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PPAIRS = PPAIRS - 1 
GO TO 36 

ENDIF 
EfIDI F 
GO TO 1 

40 CONTINUE 
QSIZE = 1 
QPART(l) = FPT 
QTIME! QSIZE) = CLOCK 

I F( QMAX! QRUN) .LT. QSIZE) QMAX! QRUN) = 1 



C. FRP-L(H) Modifications 

* READ IN FAILURE RATES 
DO 1 J 1 , TOTPRT 

READ (2,2) LAMBDA(J) 

* LAMBDA! J ) = 0.0167 

2 FORMAT ( F3 . 5) 

1 CONTINUE 

DO 3 J = 1, TOTPRT 

READ (3,8) SRATECJ) 
x SRATE(J) = 0.84 

8 FORMAT ( F8 . 5) 

3 CONTINUE 

SUMLAM = 0 

DO 7 J = 1, TOTPRT 

SUMLAM = SUMLAM + LAMBDA! J) 

7 CONTINUE 

DO 9 1=1, TOTPRT 

FACTOR(I) = PLANES X LAMBDA(I) X TOTMSN 
SPARE(I) = NINT( FACTOR! I ) + K X SQRT! FACTOR! I )) ) 
IF! SPARE! I ) .LT. 0) SPARE(I) = 0 

9 CONTINUE 



X ASSIGN PRIORITY TO PARTS VIA SUBR. QPRIOR 

CALL LINEUP(LAMBDA, TOTPRT) 

* CALL LLRAND II, GENERATE RANDOM NUMBER ARRAYS 
CALL RANDOMCA, AA, AAA) 

* INITIALIZE A NEW MISSION 



xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
X THE QUEUES 

30 QCUST = QCUST + 1 

IF(QSIZE .LT. 1) GO TO 40 
QSIZE = QSIZE + 1 
QPART ( QSIZE) = FPT 
QTIME! QSIZE) = CLOCK 

IFCQMAX(QRUN) .LT. QSIZE) QMAX(QRUN) = QSIZE 

DO 304 J = 1, TOTPRT 
DO 305 I = 1, QSIZE 

IF! QPART ( I ) .EQ. J) THEN 
LINEUP(I) = PRIOR! J ) 

ENDIF 

305 CONTINUE 

304 CONTINUE 



PAIRS = QSIZE - 1 
DONE = 1 



301 I F ( DONE .EQ. 1) THEN 

DONE = 0 

DO 302 I = 1, PAIRS 
IF(LINEUPCI) .GT. LINEUP(I+1)) THEN 
TEMP = QTIME(I) 

TEMPQ = QPART(I) 

TEMPQI = LINEUP(I) 

QTIME(I) = QTIMECI+1) 

QPART(I) = QPART ( 1 + 1 ) 

LINEUP(I) = LINEUPCI+1) 
QTIMECI+1) = TEMP 
QPART ( I + 1 ) = TEMPQ 
LINEUPC 1+1 ) = TEMPQI 
DONE = 1 

ENDIF 

302 CONTINUE 

PAIRS = PAIRS - 1 
GO TO 301 
ENDIF 

DO 303 I = 1, QSIZE 

303 CONTINUE 
GO TO 1 

AO CONTINUE 
QSIZE = 1 
QPART ( 1 ) = FPT 
QTIMEC QSIZE) = CLOCK 

IFCQMAXCQRUN) .LT. QSIZE) QMAXC QRUN) = 1 
GO TO 1 



xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



SUBROUTINE L I NEUPC QLAMBD, QTOTPT ) 

INTEGER QPRIORC 100 ) 

REAL QLAMBDC10), QTOTPT 
SAVE QPRIOR 
COMMON/ RULE/QPRI OR 

DO 1 J = l, QTOTPT 
QPRIORC J ) = J 

1 CONTINUE 

PAIRS = QTOTPT - 1 
DONE = 1 

2 I F( DONE . GE. 1)THEN 

DONE = 0 

DO 3 J=l, PAIRS 

I F ( QLAMBDC J ) . GT . QL AMBDC J + l ) ) THEN 

TEMP = QPRIOR(J) 

QPRIOR(J) = QPRIORCJ+l) 

QPRIORC J + l ) = TEMP 
DONE = 1 

ENDIF 

3 CONTINUE 

PAIRS = PAIRS - 1 

GO TO 2 

ENDIF 

RETURN 

END 



r>5 



LIST OF REFERENCES 



1 . 



2 . 



3 . 



4 . 
5. 



Naval Air Systems Command Report NAVAIR 41032A, OLSP NO. 

AV- 199 , Opera tipnal Logistic? Support Plan. P-3 

Intermediate Maintenance Advanced Base Functional 
Component (ABFc) . Revision F, February 1987. 

Rand Report R-2886-AF, The DYNA-METRjC Readiness 
Assessment Model; Mo t ivation. Capabilities, and Use , by 
Raymond Pyles, July 1984. 

Naval Postgraduate School Report 



Naval Postgraduate school Report NPS55-81- 



NPS55-81-005, The New 

' age-- 



Ross , S . M . , Introduc tion to Probabil ity Models . 3d ed.. 
Academic Press, Inc . , 1985 . 



Larson. H.J 



Introduction 

Inf eren c e., 3d 



ed . , 



£r 




and 
ons , 



66 



INITIAL DISTRIBUTION LIST 



1 . 

2 . 
3. 

4 . 

5. 

6 . 

8 . 

9 . 
10 . 

11 . 
12 . 

13 . 



No . 



Defense Technical Information Center 
Cameron Station 

Alexandria, Virginia 22304-6145 

Library, Code 0142 

Naval Postgraduate School 

Monterey, California 93943-5002 

Lt. Peter J. Latta, USN 

3062 Camino Aleta 

San Diego, California 92154 

Mr. Irving Cohen 
Karen Isaacson 
T. Lippiatt 
RAND Corporation 
1700 Main st . 

Santa Monica, California 90406-2138 

Professor D.P. Gaver, Code 55 Gv 
Department of Operations Research 
and Administrative Sciences 
Naval Postgraduate School 
Monterey, California 93943-5000 

Professor P.A. Jacobs, Code 55 Jc 
Department of Operations Research 
ana Administrative Sciences 
Naval Postgraduate School 
Monterey, California 93943-5000 

Commanding Officer 
VP-31 

Naval Air Station 

Moffett Field, California 94035-5008 

Commander 
Patrol Wing 10 
Naval Air Station 

Moffett Field, California 94035-5022 
Commander 

Naval Air Systems Command 

ATTN: PMA 240 (P-3 Program Office) 

Naval Air Systems Command Headquarters 
Washington, D.C. 20361-0001 

Lt. Stephan B. Latta, USN 
COMSEVENTHFLT Staff 

F.P.O. San Francisco California 96601-6003 

Center for Naval Analysis 
Attn: Library 
4401 Ford Ave . 

Alexandria, Virginia 22302-0268 

Col. Buck Rogers, USMC(RET) 

47-337 Lulani st. 

Kaneohe, Hawaii 96744 



Copies 

2 

2 

2 

2 

1 

1 
1 
2 

1 
1 

1 



67 



HHro 






Latta 

A comparison of six 
repair scheduling policies 
for the ‘P-3 aircraft. 




Thesis 

L28133 

c.l 



A comparison of six repair scheduling po 




3 2768 000 78456 5 

DUDLEY KNOX LIBRARY *■ 





